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TECHNICAL NOTE NO. 1878 


STUDY OF UNSTEADY FLCW DISTURBANCES OF LARGE AND SMALL 
AMPLITUDES MOVING THROUGH SUPERSONIC OR SUBSONIC 
STEADY FLOWS 
By Robert V. Hess 

SUMMARY 


A study of unsteady flow disturbances of large and small amplitudes 
moving through supersonic or subsonic steady flows is presented in three 
parts. 

In part I a point— by— point method is developed- for the calculation 
of unsteady flows through tubes with variable cross section under the 
assumption of constant flow velocity at a given cross section. The 
paper extends the work done previously by giving a detailed treatment 
of the interaction of strong shocks and large temperature contact 
discontinuities and by presenting the shock calculations and the calcu- 
lations of flows with initial entropy gradients in a form convenient 
for computation by use of computing machines. Under certain assumptions 
the formulas established may also be used for the calculation of flows 
with continuous heat addition over a large space. 

In part II calculations are made of the flow pattern created by 
the bursting into a vacuum of a diaphragm at the minimum section of a 
supersonic nozzle without a second throat. The transition time from 
the starting of the flow to the attainment of approximately steady 
flow conditions is sufficiently short to permit the use of very-short— 
duration testB. The transition time for the specific nozzle is 
presented in such a form that a "similarity rule" can be established 
concerning the transition time for nozzles of different size but of the 
same or affine shape. 

In part III integral relationships invariant with respect to time 
are developed which describe as a whole the behavior of unsteady flow 
disturbances of large and small amplitudes. The invariant integrals 
are the conservation laws for the mass, energy, potential, and pulse 
area. For the special case of disturbances of small amplitude and 
small length traveling through tubes with small cross-sectional 
gradient (short disturbances), growth and reflection of mass, 
energy, and pulse area of a disturbance traveling through a steady 
flow in variable cross section can be separated and presented as a 
function of the steady-flow Mach number. The calculations show the 
interesting result that the conditions for zero reflection of mass. 
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energy, and pulse area exist at a steady— flow Mach number 

Mo = ]Jy~ Y (which equals {5 for air), where 7 is the ratio of 

specific heats, rather than at Mq = 1 ; they also show that for practi- 
cal purposes for the range of Mach number from 1 to 5 (for air) the 
reflections are small enough so that the mass, energy, and pulse 
area of the original disturbance may be considered constant. At 
low subsonic and at high supersonic Mach numbers, however, the 
reflections may not be neglected. 


INTRODUCTION 


The study of unsteady-flow problems has recently achieved new 
significance in aeronautics, partly because of the possibilities of 
using unsteady-flow phenomena for the improvement of the performance 
of high-speed internal-flow systems and wind tunnels and partly because 
of the greater understanding of the basic nature of flow phenomena 
such a study has to offer ♦ 

A few of the unsteady-flow problems that are useful to the 
aeronautical engineer are indicated. The study of the stability of 
shocks in diffusers involves a whole series of problems. One such 
problem is concerned with the stability of a normal shock with respect 
to disturbances moving upstream; such disturbances may be produced, 
for example, by the fluctuations in a burner of a Jet airplane or by 
other types of fluctuations occurring in the operation of a Jet airplane 
or of a wind tunnel. Still another type of stability problem deals with 
the two possible equilibrium positions of a shock in front or inside of 
a diffuser; the possible "Jumping" of the shock from one equilibrium 
position to the other in steady-flow terminology is actually an unsteadv- 
flow phenomenon. 

Aside from these stability problems a series of problems involving 
flow discontinuities other than shock waves are of interest. One type 
of discontinuity is the temperature contact discontinuity created by 
the sudden pressure increase in a fluid such as occurs, for example, by 
bursting a diaphragm, detonation, or the crossing of two unequally strong 
shocks. Other types of discontinuities are the flame front moving through 
a combustible gas or the condensation front produced in a wind tunnel. 
Aside from the problems of discontinuous heat addition, there also exist 
problems of continuous heat addition over a large space (for example, 
those due to combustion) . Naturally, both cases, the discontinuous flame 
front and the continuous heat addition over a large space, are only 
idealized models of the combustion process (though often very useful) . 
Finally, for stability problems and related problems the study of unwii 
or infinitesimal disturbances moving through the subsonic or supersonic 
portion of a wind tunnel or through the internal-flow system of a Jet 
airplane is of importance. 
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The present paper is presented in three parts. Part I may he 
considered a manual in the procedure for the construction of unsteady- 
flow patterns involving strong shocks and large temperature contact 
discontinuities under the assumption of constant flow velocity at a 
given cross section (one-dimensional flow). The flow variables are 
given in simple expressions in a form convenient for computation by use 
of computing machines. The emphasis in part I is on the method of 
performing the calculations, rather than on the elucidation of all their 
physical and mathematical meanings. The method is an extension of the 
methods developed in references 1 to k. 

In part II an example is calculated for the motion of a disturbance 
of large amplitude through a tube with variable cross section. The 
large— amplitude disturbance is produced by the bursting into a vacu um 
of a diaphragm at the minimum section of a supersonic nozzle without a 
second throat. The unsteady-flow pattern thus created is calculated up 
to the time that conditions very close to steady flow are approached in 
the nozzle (the actual approach is asymptotic). The transition time 
from the starting of the flow to the attainment of approximately steady- 
flow conditions for the specific nozzle is presented in such a form that 
a "similarity rule” can be established concerning the transition time 
for nozzles of different size but of the same or affine shape. 

In part III integral relationships are developed to describe as a 
whole the behavior of unsteady flow disturbances of large and amni i 
amplitudes traveling in steady subsonic or supersonic flows through tubes 
with variable cross section. The development of such relationships is 
possible, in spite of the fact that, because of the inclined walls, a 
very complicated pattern of deformations, reflections, and re-reflections 
may occur within the disturbance, since integrals over the whole distur- 
bance exist which are independent of time. First, one such integral 
relation is derived that represents the conservation of the velocity 
potential of an isentropic disturbance. Since the expression for the 
potential contains the first power of the velocity, it is obviously 
related to the momentum, which, in contrast to the tube with constant 
cross section, is not conserved because of unknown pressures at the 
inclined walls. The other two conservation laws for mass and energy 
can be used directly. The value of such invariant integrals is doubtful 
as long as the original disturbance, its reflections, and re— reflections, 
still belong to the integrand. Therefore, as a first step for the sepa- 
rating of the original disturbances and their reflections, the concept of 
pulse area is introduced. The reason for the usefulness of the pulse 
areas is that they have as. amplitudes certain linear combinations of 
the flow variables, that is, the parameters of the families of character- 
istics, which are associated with disturbances traveling with the speed 
of sound to the right and to the left relative to the fluid. The final 
step in separating the original disturbance and its reflection is 
achieved by restricting the length and the' amplitude of the original 
disturbance and the inclination of the tube walls to emal 1 values 
(short disturbance). The investigations in part HI are an 
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extension and an elucidation of the studies made in references 5 and 6. 
In reference 5, the concept of pulse area for short disturbances is 
introduced without relating it to the potential and the momentum of 
disturbances of large and small amplitudes. The calculations are 
restricted to a short disturbance moving upstream in a subsonic diffuser 
in the proximity of a steady— flow Mach number of 1.0, and the result is 
obtained that, although the shape of the original pulse degenerates, 
its area is constant within the accuracy of the approximation used in 
reference 5* When the disappearing inclination of the walls is 
considered, this result can be misinterpreted to signify that the area 
of the short original pulse has a constant value near the minimum cross 
section of the nozzle. In part III, however, the area of the original 
pulse is shown to have finite and unique values as a function of the 
steady— flow Mach numbers with a finite derivative of the Mach number 
function at the Mach number of 1.0. Reference 6 presents the laws for 
the conservation of mass and energy (but not of the potential and the 
pulse area) for disturbances of large and small amplitudes in terms 
of the parameters of the characteristic families but is not extended 
to the case of short disturbances. The present analysis has the 
advantage of offering a clearer idea of the principal effects occurring 
in the motion of small disturbances. 


SYMBOIS 


a velocity of sound 

t time 

T temperature 

y distance along tube axis; in example of part II y is 

referred to y/Fmin 

F cross— sectional area of tube; in example of part II 

F is referred to Fm-in 

M Mach number 

v flow velocity 

Au speed of shock wave relative to flow 

u absolute speed of shock wave 

Av velocity increment through shock 

s entropy 


P 


pressure 
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n 

p 

♦ 

0 - 
7 


[i 

or 7 
P, Q__ 


dP, dQ 


m, n 

D 

d 

Subscripts: 

0 

1 

2 

o 

A,B,C 


total pressure recovery ratio 
density 

stream function 
potential 

ratio of specific heats 

parameters of characteristic families; used as 
quantities for measuring amplitudes of 
disturbances of large and Bmall amplitudes 
moving in steady flow or in gas at rest . Also 
used as mere labels for distinction between 
two groups of characteristic disturbances. 

The expressions "disturbance" and "pulse" are used 
interchangeably in present paper. 

amplitudes of growth and reflection of small- 
amplitude disturbances 

slopes 

constant 

diameter 


reference conditions 
condition ahead of shock 
condition behind shock 
steady-flow values 
at points A,B,C 


I.- A METHOD OP CALCULATING UNSTEADY FLOWS CONTAINING STRONG SHOCKS 
AND LARGE TEMPERATURE DISCONTINUITIES IN TUBES OF 
VARIABLE CROSS SECTION 


GENERAL CONSIDERATIONS 


A point— by— point method is developed for the calculation of unsteady 
flow disturbances of large amplitudes moving through tubes with variable 
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cross section under the assumption of constant flow velocity at a given 
cross section. The method used in this part is based on the. method of 
characteristics for steady rotational supersonic three-dimensional flows 
with axial symmetry developed in reference 1 and modified in refer- 
ences 2 and 3. The steady— flow method was applied to unsteady flow 
through tubes with variable cross section in reference 4. 

The present paper extends the work reported in reference b by 
giving a detailed treatment of the interaction of strong shocks and 
large temperature contact discontinuities and by presenting the shock 
calculations in a form convenient for computation by use of computing 
machines. In the course of the analysis, expressions are also given 
for the effect of the entropy gradient behind a shock created by the 
varying energy losses in the motion of the shock with varying speed 
on disturbances traveling through this entropy gradient. 

In the literature, problems of discontinuous heat addition in a 
flame front and problems of continuous heat addition over a large space 
have also been treated. The problem of discontinuous heat addition 
in the flame front (or condensation front) which has a nature s imilar 
to other discontinuity problems as, for example, the motion of a 
temperature discontinuity discussed in case C of part I has been 
treated, for example, in references 7 and 8. The problem of continuous 
heat addition over a large space has been treated in reference 9. The 
calculations in the present paper may be easily modified to include the 
case of continuous heat addition over a large space by treating the 
effects of heat addition and the compressibility effects separately. 

Since the energy equation (see, for example, equation (3) of reference 9) 
indicates that the effects of heat addition are independent of the 
compressibility effects for the case of small amounts of heat added 
at constant volume, the continuous heat addition over a large space 
may be substituted by the distribution in space of small amounts of 
heat added suddenly (at constant volume). During the time intervals 
between these sudden heat additions, the entropy along the time 
histories of the fluid particles is assumed furthermore to be constant. 
The scheme of sudden heat additions also permits a simple inclusion 
of further additions like those of fuel mass. 

A few remarks seem pertinent at this point concerning the 
applicability to practical cases of the two previously mentioned 
combustion models, that of discontinuous heat addition in a flame 
front and that of continuous heat addition over a large space. The 
first model applies to problems of combustion for which one or several 
separate flame fronts exist or where the length of the burner is small 
enough relative to the entire length of the internal flow system to 
permit the continuous heat addition being averaged by a single flame 
front without essentially changing the unsteady— flow pattern as a 
whole. The choice of proper combustion model or models for the 
analysis of the cycle of a pulse jet or of the instabilities of a 
steady— flow ram jet depends on the merits of the individual case. 
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The method of calculation given in part I is mainly "based on the 
presentation in the form of difference equations of the variation in 
convenient linear combinations of the flow variables along the time 
history of disturbances traveling with the speed of sound relative to 
the fluid (compression and expansion waves) . In mathematics, the time 
histories of the disturbances which do not exceed the order of magni tude, 
of discontinuities in the velocity gradient across the t ime histories 
are called the families of characteristics.. The characteristics have 
furthermore the significance that they also exist when no disturbances 
are moving along them. The characteristics represent thus, in effect 
the "characteristics" of the flow structure. The families of charac-* 
teristics are hereinafter called for brevity the "time histories of the 
disturbances" in contrast to "the time history of the shock" (or strong 
shock) for which the discontinuity occurs in the velocity itself. The 
linear combinations of the flow variables are called "the parameters of 
the characteristics ." They are used as quantities for measuring amplitudes 
of disturbances of large and small amplitudes moving in steady flow or in 
a gas at rest. They represent also mere labels for distinction between 
the two groups of characteristic disturbances. The labels of the charac- 
teristic disturbances refer to two different reference systems, one for 
the subsonic case and one for the supersonic case. The reason* for the 
difference in reference systems for the characteristic p aram eters in the 
subsonic and supersonic cases lies in the following: For the subsonic 

case (v < a) the fluid velocity may change between down-tube or up-tube 
direction; whereas for the sunersonic case (v > a) the fluid velocity 
may have only one direction. (The down-tube and up-tube directions or 
positive and negative y-directions may be arbitrarily designated in the 
upstream or downstream directions.) Thus, for the subsonic case the 
proper reference system for the distinction between the two groups of 
characteristic disturbances is one at rest with respect to the tube; 
whereas for the supersonic case the proper reference system is one 
moving with the fluid velocity v. In most literature on the subject 
of unsteady flow these parameters are represented, with the use of 
different signs, by the pair of symbols X and n, P and Q, or r 
and s. In part I the symbols X and p are used in agreement with 
reference 4 . As well as being identified with the motion of a small 
disturbance, the characteristic parameters have the significance that 
they remain constant along the time histories of the disturbances for 
isentropic flow through constant cross section. For the general case of 
unsteady flow through variable cross section with entropy gradient and with 
or without gradual heat addition, the time histories of the disturbances 
moving relative to the fluid (compression and expansion waves) are not the 
only characteristic families of the flow. Since the entropy is integrable 
or constant along the time histories of the fluid particles, these time 
histories form a third family of characteristics. 


METHODS OF CALCULATIONS 
General Remarks 


, , Unsteady-flow problems through a tube deal with the variation of 
the flow variables (for example., the velocity of flow v and the velocitv 
of sound a) along the length of the tube with time. In order that the 
unsteady-flow phenomena may be determined as a function of time the 


8 


NACA TN No. 1878 


flow conditions along the tube must he known for a given time (initial 
conditions). In addition, the changes in flow boundaries, including 
flow changes at ends of the tube (boundary conditions), must also be 
known as a function of time. In most problems considered, the cross 
section does not vary with time and, thus, the boundary conditions 
apply only to the ends of the tube. 

The calculations are applied to the problem of one-dimensional 
unsteady flows through tubes of variable cross section. The assumption 
of one-dimensional flow is made on the basis that the rate of change 
of cross-sectional area with respect to distance is small and that the 
flow velocity at a given cross section can be assumed constant. The 
tube cross section is assumed not to vary with time. As in previously 
developed methods, the effects of friction and heat transfer have been 
neglected. (The physical nature of these effects does not permit 
treatment with the theories dealing with wave propagation.) 

By application of the method of characteristics to the three basic 
flow equations, the equations of the characteristic families and of the 
variation of the most significant combination of the flow variables 
along them are established. For the case under consideration, the 
characteristic families consist of the two families of the possible 
t ime histories of the disturbances propagated with the speed of sound 
relative to the fluid and the single family of time histories of 
convective variations in entropy transported with the velocity of the 
fluid particles. Since the purpose of the calculations is partly to 
determine the equations of the time histories of the disturbances and 
the convective entropy variations, naturally these equations cannot 
be given in advance. The slopes of the unknown time histories, at any 
point of the flow field, are the velocities of the disturbances and of 
the convective variations in entropy (or of the fluid particles). The 
particle velocity at a given time t expressed in the convenient 
y,aot coordinate system is given by 

1L = _L (1) 

ao dt aq 

The velocity of a downstream-moving disturbance at a given time t is 
given by 

= v + a (2) 

ao dt ao 

The velocity of an upstream-moving disturbance at a given time t is 
given by 


1 dy v — a 
a 0 dt ‘ a 0 


( 3 ) 


The variation of the most significant combination of flow variables 
along the three families of characteristics is given next. For the 
single family of time histories of convective variations in entropy, 
the entropy itself is naturally the significant combination of the flow 
variables. For the families of the time histories of the disturbances 
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traveling relative to the fluid, the most significant combination of the 
flow variables is represented by the parameters \ and u, where 


\ = 







(k) 


and 


U 




( 5 ) 


The method of calculation is developed and applied to the following 
three cases: 

Case A — Unsteady flow with an entropy gradient in a tube of 
variable cross section 

Case B — Unsteady flow containing a strong shock through a tube 
of variable cross section 

Case C — Unsteady flow containing a strong shock and a large 
entropy or temperature contact discontinuity through a tube of variable 
cross section 


General Formulas 


Guderley (reference 4 ) resolves the problems of unsteady flow into 
one of integration of a system of three partial-differential equations 
by means of the method of characteristics. The three equations are: 
the equation of motion 


I|P + V |2 + |X =0 

P oy dy dt 


the continuity equation 






pv 


d log F 

ay 


ana the energy equation 
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Since Guderley only calculates the cases for which an entropy gradient 
exists "behind a shock due to varying energy losses in the shock, the 
energy equation states that the entropy is constant during the time 
history of a fluid particle. 

By application of the method of characteristics to the three basic 
flow equations, the variation of the parameters X and p along the 
time histories of disturbances moving in an entropy gradient produced 
behind a shock because of varying energy losses in the shock is obtained 
as (see references 4 and 10 ) 


and 


1_ dX = jz_ d lQ g l , _a_ ds 
a 0 dt a 0 a 0 dy 7 B dt 


1 d M- a_ _y_ d log F _a_ 4s 

a 0 dt a o a 0 dy 7 B dt 


(6a) 


(6b) 


or in terms of the stream function \|r and the total pressure 
recovery IT 


2 7 

dX _ __a_ v_ d log F _ i/a_\7-l dll 

a 0 dt a o a 0 dy 7 \ a 0/ di 


and 

\ 


27 

_1_ a_ _v_ d log F 1 / a -p 

a 0 dt a 0 a 0 dy 7 V a 0/ d\|r 


(7b) 


From equations (6) and (7), the parameters X and n can be seen 
to have the added significance that they remain constant along the time 
histories of disturbances for isentropic flow through tubes with constant 
cross section. (See also reference^.) The physical significance of 
this constancy is that if one chooses the characteristic parameters 
as the amplitudes of the disturbances, the amplitudes of disturbances 
moving downstream and upstream through tubes with constant cross section 
do not interfere with each other. * 
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The stream function i is defined at each point C of the 
diagram of y against a Q t "by the line integral 


♦ = 




1C 

uo 


(r ^ dy - r ^ t it) 


(8) 


Equations ( 6 ) and ( 7 ) differ by the factor F 0 from the 
corresponding equations of reference This difference is due to 
the fact that Guderley uses for the stream function the dimension of 
a length as is frequently done for the calculation of steady supersonic 
flows. For the present calculations, using a volume as the dimension of 
the stream function seemed more convenient; the factor Fq was thus 
eliminated. The parameters IT and are both constant along the time 
history of a particle between the shocks. The parameter U for the 
total pressure recovery of the shock is given by 


n = 


2_ 

/a 2 \ P2 

\ a l/ P 1 


( 9 ) 


The quantity p/p 0 appearing in the expression for \|r (equation ( 8 )) 
is consequently related to II by 


2 

7-1 


Po = n (i%) 

For the purpose of setting up the computational equations, the 


( 10 ) 


v a v + a v — a 


quantities , , 

ao a o a 0 a 0 a 0 a 0 dy 

may be expressed in terms of the parameters X 


_a_ _v d log F ^ 1 /js^- 1 F m 

A /■> A H V * 71 O. ) A a »!<- 


A*o) 


dt 


and p as follows: 


v _ X — p 

ao 2 


( 11 ) 
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fo = 1 + 2L i r l(x + ^ 


( 12 ) 


= 1 + 21±J: X - 3_Z_Z ^ 

a 0 4 4 


(13) 


= -1 - h' + x 

a o 4 k 


(14) 


a v d log F 
a 0 a 0 dy 


■ t 1 + ^ ♦ “>](^ I g) <«> 


and 


27 

7-1 


27 

7-1 


F ^-rh 2 r 1(l + u) ]~ F 


dH 

d\|r 


(16) 


In order to calculate the flow field, equations ( 7 a) and ( 7 b) are 
evaluated in difference form; this evaluation permits the construction 
of the unsteady flow field "by means of a point— by— point method in the 
following manner: 

The flow variables v/a q and a/aQ are known for two locations 
in the tube y A and y B at two known times a Q t A and a Q t B (see 

fig. l) . From this information, the flow variables are calculated at 
a third location in the tube Yc at a corresponding time agtQ. The 

location Yq, the corresponding time &QtQ, and the values of the flow 

variables at C are determined by the intersection of the time— history 
curves of the waves through A and B and the knowledge of the variation 
of the flow variables along these time— history curves. Since the 
curvature of the time histories, however, depends on the variation of 
the flow variables along them, an iteration process was used in 
references 4 and 10 for the simultaneous determination of y c , a 0 t c , 

and the flow variables ^c / a 0 an ^ a c/ a 0 * 

In reference 4 the known time— history curves of disturbances 
traveling through constant, cross section were used as a first approxi- 
mation. (The curvature is due to interference of the two families of 
time histories of disturbances.) Since, however, the curvature of 
time histories of the waves traveling through a tube of variable cross ■ 
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section is still unknown, the use of the curvature for constant cross 
section as a first approximation introduces unnecessary complications. 

The method of the present study is "based on the substitution of 
the very slightly curved time histories through the points A and B 
"by the tangents at these points as a first approximation. In this 
first approximation, the point C obtained by the intersection of 
the tangents is used for the calculation of the flow variables. Within 
the order of accuracy of the calculations, the point C obtained by 
intersection of the tangents is identical with the point obtained by 
the intersection of the time— history curves in the first approximation 
of reference 4. The second approximation consists in using the 
arithmetic mean between &A/ a 0 and a c/ a Oj yj^/slq and vq/bq, a B/ a 0 

and ac/aCb and Y B/ a 0 and vp/ao for determining new time histories 
between A and C* and B and C’, respectively. The values sq'/sq 

and v c '/ a o for " the second approximation are found by integrating along 
the newly found time-history curveB. 

It should be emphasized at this point that the use of a second 
approximation does not imply the UBe of very large steps in the 
construction of the net of time histories, since, for the success of 
the second iteration, terms of the second order must be small. If 
small steps are used initially, however, the use of a second approximation 
may be unimportant by engineering standards. Also, for the use of more 
elaborate computing machines such as the Bell Telephone labora- 
tories J.-66'jkk relay computer in use at the Langley Laboratory, the 
first approximation used with many very small steps may prove more 
advantageous . 

In developing the present method, the possibility of directly 
obtaining more exact values at C was also considered. Th 9 se values 
could be obtained by intersecting lines AC and BC (fig. l) given 
by the average directions between A and C and B and C, 
respectively. Furthermore, the average conditions between the flow 
variables at A and C and B and C, respectively, could also 
be taken into account. The use of average conditions, however, 
results in cubic or quartic equations for the unknowns instead of the 
linear equations. The solution of these highen-order equations would 
itself require" an iteration process. 


Case A — Unsteady Flow with Entropy Gradient in a 
Tube with Variable Cross Section 

The purpose of the following calculations is to determine the 
effect of unsteady flow disturbances produced isentropically by 
compression or expansion on the 'change of the flow variables with 
tim 9 as the disturbance travels in a tube with variable cross section 
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through the entropy gradient "behind a shock. With the aid of assumptions 
stated in the section entitled "General Considerations," the expressions 
lend themselves also to the calculation of problems of continuous heat 
addition over a large space. The calculations are presented conveniently 
in the following manner (see fig. l) : 

Two points in a tube with given variation of cross section F = f(y) 


A I a + rr 1? ^A T A ^ -i— j- . ' 

T° 7aj Aj v v \ j 11 A > 


and 


*(•<** »B. F B, "B, * B ) 


are given. The point 


'fo l P ^ V r 0 ’ r 0 ’ n 0’ *c) 


is to be determined. 

The intersection of the tangents to the time— history curves 
through points A and B results in 


a o t c 


e - d 
l — k 


and 


y c = 7 a 0 t c + d , 

where the constants in terms of X and 4 (equations (4) and (5)) 



are 
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d = -la 0 t A + j A 


e = -ka Q t B + y B 


The quantity y c determines F c and dF c /dy c since the shape of the 

tube is given. The next step is to express equations (7a) and (7b) in 
difference form. Based on the first approximation of tangency, X c 

and (j^ are given in the following form: 


l c - X A 


1 ZA 1 1 dF c \ 

2 a 0 a 0\^A ^A F c ^C/ aoA CA 


27 

1/^aV" 1 f c + f a n c ~ 
7 \ a 0 / 2 t c - t A 


a O At CA 


(17) 


and 


fiQ = hs - 


1 *g/l_ !ES 

2 a 0 a 0V F B d ^B 


+ ^ S§) aoAtcB 


27 

A 1/^B V -1 F c + f b 

7 \ a 0/ 2 


n c ~ n B 
tc ~ t B 


a O At CB 


(18) 


The flow variables v c /a 0 and ac/a 0 are determined from X c and p c 
by means of equations (4) and (5). 
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All quantities on the right-hand side of equations (17) and (18) 
are now known, except Hq and which may he obtained by assuming a 

linear variation of the values of II end \|r from A to B. Since U 
and i|r are constant along the time-history curve of a particle, and 
because the time— history curve of the particle is here substituted for 
the tangent, the condition of linear variation of II and \|r along AB 
yields 


n B ~ _ ^B ~ 

n c - H A " i|r c - * A 


fl 

c 2 


The fact that Hq and i|fQ are subtracted in the n ume rator is 

necessitated by the physical condition that the velocities of particles 
flowing through BC have to be of the same sign as the velocities of 
the particles flowing through CA. Since the tangent to the t ime -history 
curve of the particle bisects the angle formed by the tangents to the time- 
history curves of the disturbances at C, the ratio c l_/ c 2 . i 0 equivalent 
to the ratio b/a. (See fig. 1 .) The quantities IIq and i|Tq are 
then given by th 9 following formulas: 


n c = 


Uga + n A b 
a + h 


= + *A* 

^ a + b 


where 


= \/( a 0 t C - a 0 t A) 2 + (y.C - ^a) 2 


and 
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Finally, a "brief note is made concerning the scheme of sudden heat 
additions (at constant volume) distributed over a large space. The 
sudden additions of heat are made at the intersections of the net of 
time histories of the disturbances. During the time intervals between 
these sudden heat additions, the entropy of the fluid along the time 
histories of the fluid particles is assumed furthermore to be constant. 
As is shown in figure 1, the time interval between the sudden additions 
of heat is the time it takes a particle t,o move albng its time history 
from the connecting line of the two former net points to the new point. 
The variations of X and n along the time histories of the downstream 
and upstream moving disturbances are calculated by the use of 
equations (6a) and (6b) in difference form. 


Case B — Unsteady Flow Containing a Strong Shock through 
a Tube with Variable Cross Section 

The purpose of the calculations in this section is to determine the 
effect of a strong shock on the change of the flow variables with time 
as the shock travels in a tube with variable cross section through a 
steady isentropic flow or through a gas at rest. The more complicated 
case of a shock crossing an unsteady isentropic flow is not developed 
herein. A method very similar to that used for the "interweaving" into 
the flow of a temperature contact discontinuity (see case C) can, 
however, be used in that case. 

The expressions for the flow variables are developed for the case 
of a shock moving upstream through a steady isentropic flow in a tube 
with variable cross section. The presentation of the flow problem is 
given in figure 2. A steady flow is assumed to exist in the direction 
of the negative y— axis and a shock is assumed to be produced traveling 
in the direction of the postive y— axis (through decreasing cross section) 
The production of a shock of given strength is identical with the 
condition that the following values are known: Au A ^a A ^, Av^ya^, 

and &A 2 ^ a A;L* These values in turn are used to determine the absolute 

speed of the shock u A ^aQ and the conditions behind the shock v^^a 0 

and a.p^ja.Q by the following fundamental relationships: 

u A v Aq | Au a 

a 0 a 0 a 0 

where v A /sq is negative and Au A /aQ is positive 


v A 2 _ V A X Av a 
a 0 a 0 + a 0 
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where Avj^J&q is positive, and, finally, 

a Ag _ a Ag a A]_ 

a 0 a Aj a 0 

The quantities Aua/sq and Av^/ao obtained by multiplying Au A /aA.]_ 

and Avj^/aA-^, respectively, by a A B / a O 

In order to find the changes in the behavior of the shock as it 
travels in the positive y— direction, a tangent to the time— history 
curve at B is drawn and brought to an intersection with u^/ag at 

point C. The position of the point B in the diagram of y 

against agt is known, as well as the flow variables v B /aQ and a B /a q 

at B. The point B may, generally, be assumed to lie on the time- 
history curve of a disturbance of previous construction. In this 
particular case, which presents the initial calculations for the shock 
movement in a diffuser, the point B is assumed to lie on the t ime — 
history curve of the particle v/^aQ. 

The calculations of the flow variables at the point C are based 
on the knowledge of their variation during the time history AC of the 
shock and the time history BC of the disturbance. The variation of 
the flow variables during the time history BC of the disturbance is 
given by equation (17) since the small disturbance is moving down tube, 
the flow is subsonic, and the parameter X is chosen to be associated 
with the down— tube direction. The variation of the flow variables 
during the time history of the shock is developed as follows: 

Two points 


A 



v F ^1 lAg fAi ^2 ^ 

Ja ’ a ’ dy A ' a 0 ' a 0 J a Q * a Q J a Q > 


IT 


A> 



and 


B 



f b> 


AFb V B a B 
dy B ’ a 0 ' a 0 * 





are given. The point 


a oV 




c> 


dP 0 

dye" 


I2i 

a 0 ' 


vp ap ap 
'“ / 2 ^2 

*0~’ ao”' 


^C 

& 0* 


XI C J *c 
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is to "be determined. In the present method, the location of point C in 
the aot,y plane is given by 


a otc 


f - e 
h - d 


and 


y c - ta Q t c + d 


where 


h 




= 1 + 


7 + 1 

k 



^B 


e - — ' ba Q t A + y A 


and 


f = -da 0 t B + y B 

. The location of point C, given hy y c , determines F c and dF c /dy c 
from the shape of the tube; furthermore, if Jq is known, the quanti- 
ties V C-^J ao 811(1 a C 1 / a 0 are determined since a given steady flow 
exists in the, tube. Only the determination of 
and uq/sq thus remains. 


v c 2 / a o> & c 2 l &0 ’ 
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In order to relate the values of the flow variables at A and C, 
a relation is set up between the shock parameters Av/aq and a 2 /aq; 
the relation is obtained with the aid of the shock equations (refer- 
ence ll). Thus, 


Av - 2 /Au _ a l\ 

aq 7 + l(aq SI J 


(19) 


and 


f2 

a l 


± ( 2(7 - 1 ) 

~ x 2 "1 

f—\ — 1 

/ a l\ 2 
7 + 1 — | 

(7 + l) 2 

[V a i/ 

\Au/ 


( 20 ) 


A plot of against Av/aq is given for air (7 = 1.4) in 

figure 3 and more extensive values are presented in table I. 

For the determination of the change of the flow variables along 
the time— history curve of the shock, the variation of the shock 
parameters Av/a^ and Aa/a (where Aa = — a-jJ is assumed to be 

linear. Thus, 


_ £fA _ / Av C Av a\ 

a C x a A 1 ^C-L a A J 


( 21 ) 


Since 

Aa^ = aAg _ 
a Aq a Aq 


and 

Aac _ aCg 

a Cq a Cq 


equation (2l) may be written as 





m 


( 22 ) 
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The change of the flow variables during the time history of the shock 
may be conveniently expressed in the form 


a Co v Co 

ao = m A ag + Constant 

This expression is obtained by multiplying equation (22) by aQ^jaQ. The 
following equation results: 


a C 2 V C 2 

— = m A — - 


a 0 


a 0 


R 


(23) 


where 


▼Ci a C^ 


R = m A ^ m A — a & 

& r\ n ^ \ ^ « a o a 


a 0 


a 0 


Avi 


^Ap 


a A n a A. 


For the special case of a shock moving through variable crosB section 


with air at rest in front of it 
equation (23) becomes 


( aA i - a <=i - 


V Aq V Ci 

ao and — — = — -=4 = 0 
u a 0 a 0 ’’ 


a C 2 

a 0 


v c 2 

a 0 


= m A - R 1 


(24) 


where 


*1 


m A 


7 A 2 _ a Ag> 
a 0 a 0 


The linearization in the shock equations is equivalent to substituting 
for the curve in figure 3 a tangent with the direction m for a point 
with the coordinates Av/a^ and ap/aj. . For the purpose of determining 

accurate values for m, an analytic expression for m which was 
developed from equations (19) and (20) is given in appendix A. 

A plot of m against A.v/a q for air (7 = 1.4) is presented in 
figure 4 with more extensive values given in table II. The error 
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involved in the linearization can he seen from table II to he small for 
small steps. The lower limit of m is identical with the value of m 
obtained from the laws for isentropic disturbances 


f2 7-1 Av . -1 

a l 2 a 2 + • L 



for air (7 = 1.4) 



whereas the upper limit is reached asymtotically as Av/ay approaches 

infinity. The upper limit given in reference 11 is \j ~ ~ 2 1) (0,52915 for 

air). This upper limit, however, is reached asymtotically only for a few 
practical cases, since 7 may vary greatly across the shock. 


The variation of the flow variables during the time history BC of 
the disturbance can be expressed in terms of the same two unknowns vq j&Q 

and aQ 2 ^aQ. From equation (17) and equation (4), the following expression 
is obtained: 



2 fB _ 32 

7 - 1 a 0 a 0 


1 fB 2l/l_ 

2 a 0 a 0 \F B dy B + 


1 ^c\ 

F C dy c / aoAtcB 


1 f B + F C 
7 2 



n c - n B 
ic ~ ^B 


a O At CB ( 2 5) 
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Finally, 
quantities, 
in the form 


and IT C have to he expressed in terms of known 
For the purpose of expressing \|r c , - ijr-g i s written 

~ ^A + (^A ~ ^b). In accordance with equation (8) 


*C - *A = 


F* + F P PAi + PC 


2Pq 


* 


■ ■ 


2a, 


0 


a cAtAC 




(26) 


and f A and are known from previous calculations. For the special 

case of a shock moving into a gas at rest 


, . Ffl + Fp 

+c - +A = — - g C A^AC 

The only value which still must he expressed in terms of known 
quantities is He. If a linear relation is assumed, therefore, 
between Ha and ap^jeLA^, and Hp and ac^ac^, then 


nc = "(S ■ *£) + n * 


where n is the slope of the curve H = f(a 2 /ai) for 

and H a . The derivation of the formula for the slope n is given 
in appendix B. A plot of H = f^/ai) is given in figure 5 and the 

values are given in table I. A plot of the slope - = n is 

d (a 2 /a 1 ) 

given in figure 6, and the variation of n with a 2 /ai for a shock 
is given more completely in table III. 



Substituting equations (27) and (23) into equation (25) gives finally 


2 ^ 
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and 


V °2 _ 1 a C2 R 

a 0 m A a 0 m A 


Case C — Unsteady Flow Containing a Strong Shock a Large 
Entropy or Temperature Contact Discontinuity 
through a Tube with Variable Cross Section 

The purpose of the calculations is to determine the interaction of 
a strong shock and a large entropy contact discontinuity or temperature 
contact discontinuity as they travel through a tute with variahle cross 
section. Such temperature contact discontinuities are, for example, 
produced by bursting of a diaphragm in a tube. The bursting causes* 
a shock and "centered" expansion disturbances, due to instantaneous 
expansion, to travel in the tube. The strength of the shock and the 
extent of the expansion are determined by the condition that equal 
pressures and velocities establish themselves in the flow between 
the shock and the expansion. A temperature contact discontinuity 
occurs because the expansion lowers the temperature; whereas the shock 
raises the temperature and, thus, causes two flow fields of different 
temperatures to be in contact. In figure 7 the results of bursting a 
diaphragm are presented in a diagram of y against aot. 

In the calculations presented for case A, relations had to be 
established for the variation of the flow variables along the time- 
history curves of the disturbances moving upstream and downstream. In the 
calculation of the temperature discontinuity, the variation of the flow 
variables along the time history of the discontinuity must also be 
known. For this purpose the condition that the entropy is constant 
during the time history of a fluid particle must be considered. 

As a result of this condition, the ratio of the' sound velocities 
on both sides (l and 1»; see fig. 7) of the temperature contact 
discontinuity is constant along the time— history curve. The variation 
of the flow variables along the time— history curve of the temperature 
contact discontinuity is, thus, given by: 



( 29 ) 


2 6 


NACA TN Ho. 1878 


and 


ai'/ap 

a l/ a 0 


(30) 


where D is a constant. The ratio - a / a fl is determined hy the shock 

a l/ a 0 

strength and the extent of the centered expansion disturbances . The 
relations may also be conveniently expressed in terms of X. and u 
(equations (4) and (5)).' 



(31) 


1 + * ei’) 

1 + ^ir^i + kl) 


(32) 


The flow variables along the time— history curve of the temperature 
contact discontinuity can now be calculated. (See fig. 8 in which the 
shock and the centered expansion waves at t = 0 are shown by dashed 
lines.) The following points are known from previous construction: 

Point 1 f a 0 t, , y,, F . } } 3 } il-i } ^1^ is situated somewhere 

\ dyj a o a o x V 

along the time history of the temperature discontinuity (also indicated 


in fig. 7). Point 2 (a Q t 2 , y 2 , F 2 , 11^, is situated either 

on the shock (see fig. 8) or along a disturbance traveling in the direction 

( ' dF 7 vv a 7 \ 

a 0 t 7 , y 7 , F _L — - U 7 i|r 7 1 
ill d jrj a Q a 0 ' 


dFo vq a q 

and 9 (a 0 t 9 , J 9 , F 9 , — , — , Uj, * 9 


( a 


dy 9 " a 0 * a 0 J 

history curve of a disturbance traveling in the direction of the negative 


are situated on the t iro** — 


y— axis . 


Points k (aot k , j h , P 4 , %■, n h , 


/ dF v a= 

and 5 (a Q t 5 , y^ F ? , — , 

time-history curve of the temperature discontinuity } are to be calculated. 
The values of aQt, y, and F for point 4 are identical to those of 

point 5. 


n B 


lying on either side of the 
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The first operation in the calculation is concerned with the inter- 
section of the tangent of the time— history curve of the disturbance at 
point 2 (traveling in the negative y-direction) with the tangent of the 
time-history curve of the temperature, discontinuity at point 1. The 
intersection results in 


a (A = 


and 


7k = paoV + r 


where 


P 


= X 1 ~ ^1 

2 




7 + 1 


«2 


3 ~ 7 


r = pa Q t + j 1 


and 


s = — qa^^ + y 2 


The values of and y^ are identical with the values of agt^ 

and y^, respectively. 

Use iB then made of the fact that the tangent to the time-history 
curve of a disturbance moving in the direction of the positive y— axis must 
satisfy the following conditions: 


dF 


The quantities a 0 t^ y-^, F^, — 


and must assume values 

such that the tangent to the time-history curve of the disturbance in the 
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positive y-direction at point 3 will pass through the previously determined 
point 5 j thus. 



= 1 + 


7 + 1 

~T~ 



HZll 

aQt^ — aQt^ 


(33) 


Equation ( 33 ) contains four unknowns. Three additional equations for the 
quantities aQt-^, y^, X^, and ^3 are obtained "by assuming linear 

relationships "between the quantities j, X, n, and aQt at the 
points 7, 9, and 3. Thus, 


where 


and 


y^ — maQt^ + "b 
X^ = kaQt^ + d 
= Za^yt^ + e 


(34) 


= y 9 Zll 

a 0 t 9 “ a o" t 7 


k = ^ ~ x 7 

a 0 t 9 ~ a o^7 

l = ^ ~ ^7 

a 0 t 9 ~ a o^7 

b = -ma 0 ty + yj 

(1 = “’"kO.Q'fcy + \rj 


© — “ZaQty + fly 
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Substituting the values for y 3 , X^ 3 and p 3 from equations ( 34 ) into 
equation (33) yields the following expression for aQt^ 


„ + _ -H ± \|H 2 - kGK 

a 0 t 3 = 2G 


where 


G = -2-HZ 1 — Z_t_i 


H = - 


■1 + - + 1 kantc; - 2L+2 d - i-~ 7 


3-7 


— «o^5 ~ “x“ a “ "nr* ZaQt 5 + xx 6 + m 


and 


K = 


a 0 t 5 + 


7 + 1 


a o t 5 d - y 5 


+ b - 


3-7 


a 0 t 5 e 


The quantities X X^, and ^ can now be obtained with the aid of 

the relations Jor the changes of the flow variables along the time history 
of the waves 3,5, 2 , 4 , and the temperature contact discontinuity 1 , 4 . 

+v, Th j variations of the flow variables along the time-history curves 
of the downstream and upstream waves are given by equations (17) and fl8) 
Equations (31) and (32), that is, the variation along the time-history 
curve of the temperature contact discontinuity, are written more 
conveniently as 


X 5 7 ^5 = X 4 “ ^4 

1 + 2L 1T : ( X 5 + ^)' = D [ 1 + + W*)] 


30 


MCA TN No. 1878 


The following equations are thus obtained for the two remaining unknowns: 


D + 1 X 5 D + 1(7 - 1 + 

^5 = 2 (7 — 1 + + D 2 ^ (35b) 

The quantity TT5 is on the expansion side of the temperature contact 
discontinuity and/ hence, is equal to one; in this case = 0. The 
quantity is constant during the time history of the temperature 

contact discontinuity (particle) and, hence, has the value corresponding 
to the shock which was initially produced by the bursting of the diaphragm; 
in this case \|r^ = 0. 


II.- STUDY OF THE FLCW PRODUCED BY THE BURSTING INTO A 
VACUUM OF A DIAPHRAGM AT THE MINIMUM CROSS 
SECTION OF A SUPERSONIC NOZZLE 

GENERAL CONSIDERATIONS 


As an introduction to the case of the bursting of a diaphragm in 
a nozzle, that is, a tube with variable cross section, the nature of 
the simpler problem of bursting a diaphragm in a tube with constant 
cross section is briefly discussed. A thorough account of this problem is 
given in reference 12; a statement is also made therein that seems a 
good introduction to the problem of bursting: 

On first thought, one might be led to believe that at the instanoe 
of bursting the diaphragm the total pressure Jump across the diaphragm 
would be propagated toward the low-pressure side as shock wave. However, 
this is not possible, for then the entire air mass on the high-pressure 
side of the diaphragm would have to be suddenly accelerated to the speed 
of the mass of air behind the shock. 

The gradual acceleration of the mass of air on the high-pressure 
side is accomplished by the fact that the total— pressure Jump across 
the diaphragm produces not only a shock traveling into the low-pressure 
side but also a sudden expansion spreading gradually into the high- 
pressure side. 
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In setting up the calculations for the present investigation of a 
nozzle for high supersonic Mach numbers, an infinite pressure ratio 
across the diaphragm was chosen. This is equivalent to having a vacuum 
on the low-pressure side of the diaphragm. Furthermore, the assumption 
was made that the high-pressure tank is sufficiently large to keep the 
pressure constant in it. The choice of a vacuum appears to he of 
considerable advantage since conditions for the existence of a shock are 
not fulfilled for that case, and thus, only the effect of the spreading 
of an isentropic sudden expansion of large amplitude into a tube with 
variable cross section must be calculated. 

Figure 9 gives a presentation in the plot of y (distance along the 
tube) against agt of the sudden expansion into a vacuum for a tube with 

constant cross section. The plot of y against aot is a modification 
of the length— against— time diagram which is convenient because it 
presents the velocities of the disturbances and the fluid elements in non- 
dimensional form. The lines radiating from the coordinate origin 

( y ~ aQ ~ ^ lines ) in figure 9 and the curved line, substituting for a group 

of parallel lines, ~ lines^ represent the two groups of character- 

istics. Since initially the sudden expansion has the na ture of a flow 
through constant cross section, only the radiating characteristics are 
carriers of disturbances . The group of curved parallel characteristics 
are not carriers of disturbances; however, they are necessary for the 
analysis of the flow structure. Since this problem is one of unsteady 
flow through constant cross section, one of the parameters X and n 
of the characteristic families has to be constant. In this particular 
case X will be constant since it is the down— tube characteristic which 
cuts across the infinity of radiating characteristics carrying the 
disturbances which build up the finite-amplitude disturbance. Thus, 

X indicates the variation of the flow variables in the finite -amplitude 
disturbance. If the expansion starts from air at rest and the increase 
of fluid velocities in the plot of y against a Q t is assumed to be 

positive, the slope of the time histories of fluid particles increases 
as indicated by the v/ao line in figure 9 > furthermore, for the case 
of the expansion starting from air at rest X is equal to zero through- 
out the expansion. The flow variables are then related as follows (see 
equation (4))': 




= 0 


or 


_a_ _ 1 _ 7-1 v 

a 0 2 &o 


( 36 ) 
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The speed of the disturbances traveling along the radiating 
characteristics, and therewith the slope of these time histories, 
is given hy 3 


1 d-y _ v — a 
a 0 dt ao 


( 37 ) 


It may he seen from equations (36) and (37) that the slope of the 
radiating time histories of disturhances varies from —1 to — - 

'V 


7-1 


(which equals 5 for air); for v = a the slope of the time history is 
zero. (Fig. 9, representing the initial flow conditions, is drawn in 
the scale of 2; 5 in agreement with the final steady-flow pattern, 
fig. 13 ; the final flow pattern had to he drawn in that scale to^avoid 
crowding the characteristic lines.) Aside from giving the limits 
indicated, equations (36) and (37) show that, in the subsonic domain 
of the expansion, the disturhances travel "up-tube" and in the supersonic 
domain they travel "down— tube"; whereas relative to the fluid they always 
travel upstream. The difference in behavior between the subsonic and 
the supersonic domain of the unsteady expansion can also he brought out 
well hy observing its variation with time for two cross sections, one 
located at a negative y— value and the other at a positive y— value. At 
the negative y-value (subsonic domain) the expansion will cause a/ao 
(and thus the pressure) to decrease; whereas at the positive y— values 
(supersonic domain) the expansion will cause the pressure to. increase. 


The group of 


v + a 
ao 


characteristics carries zero disturbances for the 


caBe of a single large-amplitude expansion through constant cross 
section. (The same is true for the analogous problem of a single 
expansion — Prandtl— Meyer — for two-dimensional steady supersonic flow.) 
For the case of flow through variable cross section, however, 

y + a 

the . lines represent time histories of disturbances participating 

in the flow development due to the inclination of tube walls; these 
disturbances are often called reflected disturbances. 


Since it is desirable to anticipate at least some of the results 
a discussion of the physical nature of the problem precedes -the detailed 
investigation of the calculations for the unsteady flow through variable 
cross section (nozzle). So far in this paper the behavior of the 
disturbance of large amplitude has been discussed from a consideration 
of the microscopic elements; that is, the behavior of small (characteristic) 
disturbances (the elements of a large-amplitude disturbance) moving along 
the characteristic families has been under scrutiny'. While the combination 
of these microscopic elements is simple for the case of a large-amplitude 
disturbance traveling through constant cross section and, thus, permits 
the gaining of a thorough understanding of the physical aspects for. this 
case, for the motion through variable cross section not much insight 
can be gained without making the actual calculations. The reason lies 
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in the fact that for the case of variable cross section within the large- 
amplitude disturbance as it moves past the inclined tube walls, many 
local reflections and re— reflections are created which interfere with 
each other as they travel along different families of characteristics. 

The determination of the deformation of the large— amplitude disturbance 
consequently requires a point— by— point integration process and an 
iteration procedure. For the case of constant cross section, however, 
the correlation between the microscopic and macroscopic behavior will 
be simple since no reflections are created within the disturbance. 

It is thus necessary to look to viewpoints other than the microscopic 
in order to anticipate the results of the calculation. 

A lead on some other physical aspect of the motion of an unsteady 
large disturbance may be easily obtained from consideration of the 
fact that, for the case of bursting a diaphragm in a supersonic nozzle, 
the central problem is concerned with the process of reaching steady 
flow. The problem of balance between unsteady and steady flow energy 
is the basic problem for the motion of large unsteady disturbances in 
general and not Just for the special case discussed herein. The 
anticipated flow pattern in this special case is thus the steady— flow 
pattern given by Bernoulli’s equation for one— dimensional flow. 

In order to find the proper basis for comparison with the pattern 
of characteristics (time histories) of figure 9, representing the 
initial conditions, the final steady flow has also to be interpreted 
from consideration of the time histories of small (characteristic) 
disturbances moving through it in the downstream or in the upstream 
direction with the speeds v + a end v — a, respectively. The non- 
dimensional form of these speeds is easily obtained from Bernoulli’s 
equation. For ideal gases Bernoulli’s equation can be written in the 
form: 


fo 

a 0 



( 38 ) 


or 


fa _ 

a 0 


t 


1 + 2 -r 1 % 2 


( 39 ) 
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and 


Jo 

a 0 


Mo 



m 


then 


v o * a o 


% - 1 


a 0 


i + Z -^% 2 


(M) 


and 


v o + a o 


Mp ± 1 


a 0 


f^¥ 




(*a) 


For the nozzle shape used in this paper and for air as a medium, 
the time histories of small disturbances, used to correlate the 
steady flow with the unsteady flow conditions, are given in figure 10. 

In figure 10 only three time histories are drawn, two of these are the 
time histories of small disturbances moving upstream with the speed v Q — a ( 

in the subsonic and supersonic steady flow. The third is representative of 
the time histories of small disturbances moving, downstream with the 
speed v Q + a Q . 

For the anticipation of the results of the calculations it is further 
useful to present the transition from unsteady to steady flow in a plot 
which simply demonstrates the deviations of the flow variables from the 
steady and the unsteady state. Such a plot is one with the coordi- 
nates a/a 0 and v / aQ (fig. ll) . In this plot the sudden unsteady 

expansion into the vacuum for a tube with constant cross section (initial 
condition) is represented by a straight line given by equation (36) 


_a_ = 1 _ 7-1 v_ 
a 0 2 a 0 
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Bernoulli’s equation for steady flow is represented "by an ellipse with 
the equation 



It is important for the further discussion of transition from unsteady 

to steady flow to determine the limiting conditions for — v ~ a 

a 0* a ° 

and — . From equation (36) the following values are obtained: 
maximum velocity for unsteady flow through constant cross section. 


_v_ __ 2 
a 0 . 7-1 


which equals 5 for air, and maximum velocity for steady flow 



which equals \[5 for air. Since for the maximum velocity the velocity 
of sound a is zero, the maximum of the velocity is the a**™* as that 

of — Jz ~ and — H ~ - • Further values of importance are those of the 
0 0 

critical velocities. For unsteady flow through constant cross section, 
the critical velocity is ’ 


v _ 2 

a 0 7 + 1 

which equals O.833 for air and for steady flow 


_o _ 
a 0 


\|7 + 1 


which equals 0.9129- for air. 
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The value of - corresponding to the critical velocity is zero 

in both cases. Thus in the plot of y against a Q t the "critical” 
expansion wave coincides with the agt axis. 


Two useful properties of the presentation of the transition problem 
in the coordinate system of a/ao against v/ao or the related 

system \j 2 i against v are now discussed. One of these useful 
properties is that the Mach numbers are represented by straight lines 
radiating from the center of the coordinate system. The other useful 
property concerns the geometric presentations of the speed of the waves 


for arbitrary gases as intersects of the — axis. In one presentation 

(fig. 11 ) the intersects are obtained by subtending from the abscissa v/a 

of the point A ^ the length of the ordinate a/aQ in the 

positive direction or in the negative direction depending on whether 

one desires to obtain the values of - v - p + a or that of — ~ a at the 

a 0 a 0 

point A. In the other presentation (fig. 12) the speeds of the waves 
for arbitrary gases are obtained by presenting the steady-flow ellipse 
with the equation 


2 i + v 2 = 2 i 0 ( 44 ) 

where i is the enthalpy, as a circle with the radius \j21 Q in the 

coordinate system of \j2i against v. The intersects along the v-axis 
are made for this presentation by the normals to the character- 
istics X = ^ and- H ~ • I ^ ie latter presentation 

was first given by Busemann (reference 13 ); it can also be extended 
to the presentation of the speed of shock waves in this coordinate 
system. The exact expressions for X and |i for arbitrary gases 
are given in reference 4. 

The reason for the simple geometrical picture of the speed of the 
waves lies in the following facts: The characteristic theory for 

arbitrary as well as for ideal gases is based on the fact that in the 
vicinity of each point of the flow field, the flow may be linearized 
(reference 14). Furthermore, use is made of the fact that the motion 
of the wave may be expressed in terms of a coordinate system moving 
with the wave. In such a coordinate system, the wave will stand still, 
but the fluid will be moving with a velocity v — a or v + a. Since 
the steady flow is expressed by a circle in this diagram, the variation of 
the velocity of the locally superposed steady flow is given by the motion 
of the center of the circle along the v-axis. 
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METHOD OF CALCULATIONS 
Design of Nozzle 


The shape of the nozzle is given by a simple analytic expression 
that approximates the shape of a conventional wind tunnel designed to 
give a steady— flow Mach number of 5 in the test section. With the area 
ratio thus given as 25, the length and the shape of the nozzle still 
had to be determined. The length was determined such as to give the 
shortest nozzle for the given Jfech number and was calculated to 

be 16 . 5887 ^^ n or 16.5887 & 1116 shape of the nozzle was 

given by two analytic expressions. For lack of a more exact criterion, 
the subsonic part was chosen on the basis of the reasonable criterion 
that the maximum gradient of cross section of the nozzle presented as 
a one-dimensional flow should not be larger than the gradient of 
surface area of half a spherical wave traveling into a gas at rest. 
Thus, if F and y are made nondimens ional by referring them 

to F min 8114 'ft'mln’ respectively, the variation in cross-sectional 
area is given by the parabola 


F = 2ity^ + 1 


The supersonic part of the nozzle was chosen to have the variation of 
cross section of a parabola tangent to the cross-sectional area of the 
test section and to the parabola F = 2ity 2 + 1 continued to positive 
y— values. The equation thus obtained was 


F = -Q.cms^ty 2 - 33.17786342y + 275 . 192655279 ) + 25 


The point of tangency of the two parabolas turned out to be 


F = I.33312388 


and 


y = 0.230257055 


The high number of decimal places has its basis in the fact t ha t 
a smooth Junction between the two analytic expressions requires 
considerable accuracy. This increased accuracy offers no difficulty 
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to the Bell computer since it uses a high number of decimal places 
in all computations. 


Flow Calculations 


As previously mentioned in the section "General Considerations," 
the calculations Eire concerned with the problem of transition from 
the initial purely unsteady flow condition (fig. 9) to the final steady 
flow condition. The initial conditions for the construction of the net 

of characteristics in terms of the quantities v/sq, a/ao, P, 
v a ^— , v a ~ a , and 3 = tan -1 v a ~ a are given in table IV ; their 

magnitudes Eire obtained from equation (36). In the last columns 

of the table. the values for a or tan 3 multiplied by the 

a 0 

scale factor of 2/5 and the corresponding angle tan~-*-0 tan 3^ 

are given. As previously stated, these values had to be used in the 
construction of figure 9 in order to show a basis of comparison with 
the net of characteristics in figure 13 . -The calculations were 

started at the — — a time history of the characteristic disturbance 
a 0 


which moves into bhe gas at rest 



by subtending from 


it — characteristics at an interval of O.ly. Starting with 

, a 0 . . 

these initial conditions the flow field was calculated according to 
the step—' by— step process given in part I with the aid of the Bell 
computer . 


DISCUSSION OF RESULTS OF CALCULATIONS 


The results of the calculations are presented in the plot 
of y against agt (fig. 13)., &s well as in the plot of a/aQ 

against v/sq (fig. 14 ) . In the plot of y against aQt (fig. 13) 

th 9 criterion for the attainment of steady flow is given by the 

condition that the — and ^ + a lines have to become parallel 

a 0 a 0 

to those shown in figure 10 for steady flow. It can be seen that 
this will be true for a range of higher values of aQt. The relative 

position of the two — — a lines and the -— — — line of figure 10 

a 0 a 0 

has no significance. The proper comparison between figures 10 
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and 13 is obtained by comparing the y -- ~ a lines and the — + a line 

a 0 a 0 

of figure 10 individually with th9 corresponding lines in figure 13 
by shifting them individually along the a 0 t axis. In the diagram 
of a/ao against v/ao (fig. 14) the transition from unsteady to 
steady flow conditions is conveniently expressed by the variation of 
the flow variables v/ao and a/ao at various locations y at 
the nozzle. Steady flow conditions are obtained when the y = Constant 
lines in the diagram of a/ao against v/ao reach the steady-flow 
ellipse or, in terms of Nfach number, when the Mach numbers along 
the y = Constant lines reach the steady— flow Mach number corresponding 
to the nozzle cross section at y. Since the plot in figure 14 does 
not permit the reading of values with sufficient accuracy, the values 
of v/a q, a/a 0 , and M for several values of y are given in 

table V. The corresponding values of a Q t are also given in table V 

to indicate the rate of change of the flow variables. The full 
transition from unsteady to steady flow is represented in the subsonic 
range of the plot of &./&.Q against v /sq by y = Constant lines 

starting at the point = 1, ~ = oj } since, as may be seen from 

the plot of y against aot, y, = Constant lines first intersect 

' the disturbance — = -1. In the supersonic range, the y = Constant 

lines will start at the point (3- = 5, ^ = oV since the y = Constant 

V a ° a 0 / 

lines will first intersect the disturbances = L ± a = 5 For the 

ao ao 

purpose of identification of the intersections of the -~ + £ lines 

a 0 

and the y g ~ a lines in figure 13, the lines are denoted in the 

following manner: The — a lines carry numbers ranging from —0.1 

to -2.k; the numbers refer to their initial conditions on the 

dd ne ^izlq ~ = ■^ ie V &q a ^-* nea carry angles p varying from —45° 

to 78°4l’; the angles p represent the starting angles of these lines 
for the true-ecale initial conditions. (See also table IV.) The 
scarcity of calculated points (see table V) near the maximum values in 
the supersonic region has the following reason: In the region where 

the maximum speed of the fluid = 5 ia obtained, the maximum speeds 

of the small (characteristic) disturbances — — a = 5 and v + a = 5 

a 0 a 0 ' 

will also be obtained. They are indicated by the line denoted p = 78°4l' 
in figure 13. (See also table 17.) Since the method of calculation is based 

on the intersection of time histories of — d ._ a and Y - & disturbances 

a 0 a 0 5 - 
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the size of the step of the step— by— step process will become excessive in 
the region near the maximum speed unless a very fine net of characteristics 
is used in this region. This fine net will also he necessary if higher 
iterations than those in the present calculations are used. The use 
of a not— too-fine net (used in the present calculations) will cause 

the Y + a lines and the Y - ~ a lines to become parallel before the 
a o a 0 

physical limit for parallel! ty (tan -1 5 = T8 0 41’) is reached, namely, 

at tan -1 2.3558 = 67°. No efforts were made, however, to increase 
the fineness of the net in this region since for the phenomenon to 
be investigated it was of no particular interest; and in spite of the 
use of the Bell computer, th9 use of an extremely fine net in every 
part of the flow field would cause the time required for the calculations 
to increase by factors up to 20. (The calculations of one point required 
about 6 minutes.) It should also be mentioned here that only a small 
number of the calculated points are indicated in the diagram of y 
against a Q t (fig. 13 ). 

It appears that the statement concerning the fact that the very- 
high-speed region was of no particular interest could bear a more 
detailed explanation. Comparison of figure 13 with figure 10, the final 
steady— flow picture, shows that most of the supersonic expansion 
disturbances will take no part in forming the steady— flow picture; 
they will just disappear and the steady— flow picture will be formed 
without their help. The question arises now concerning the behavior of 

the reflected ^ + a disturbances. Figure 13 shows that the reflected 
a 0 

disturbances in the supersonic region move even at a greater speed in 
the same direction (+y) as the expansion disturbances . This behavior of 
the expansion disturbances and their reflections is based on the previously 
discussed fact that the waves are produced in a coordinate system moving 
with the speed of the fluid. When the speed of the fluid is supersonic, 

both the — ~ a and the a disturbances move downstream and 

a 0 a 0 

the v - + - a disturbances have greater speed. 
a 0 


A closer investigation of figure 13 and table V also shows that 
the y = Constant lines do not completely reach the required values 
on the steady— flow ellipse. The reason lies in the following facts: 

The nozzle chosen for the calculations was a typical supersonic nozzle 
designed for parallel flow at the test section. In such nozzles the 
large variations in cross section have to occur near the mini mum 
section; thus, large changes in Lfech number also occur in this region. 

This again means that in the region near the min imum cross section large 

changes in the speed of the Y -- ~ ~ and the --■■ t- a disturbances also 

a 0 a 0 


occur. In consequence, the intersections of the 


v — a 


a 0 


and 


v + a 


a 0 


lines 
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yield inaccurate results unless very small steps are taken. The fact 
that steps taken too large could have such an effect was first noted 


during the calculations in the results obtained for the - v - + -- a time 

a 0 


histories. The excessively large steps caused the 
those with' the starting point of y = — 1 »nri 


v + a 


a 0 
aot = 1 


lines beyond 
+ a 




to become parallel and even to intersect. As pointed out previously 
in the discussion of figure 10, parallelity of the time histories 
of the disturbances would indicate that steady flow had been reached. 
It was found, however, that the reaching of parallelity was premature 
since the values of v/sq and a/aQ attained in this region were 

still far from the steady— flow values for the given nozzle shape. 


A much smaller division of the 3 lines was then taken from 

a 0 

this region on until similar inaccuracies were noted again for 

the v - + a disturbance starting at y = -2.4 / ^ - + a line -2.4V Since 

a 0 V a 0 J 

steady flow conditions corresponding to the shape of the nozzle 
almost been reached by then, . a further reduction in step size was not 
attempted. In this slightly premature stoppage of the calculations 
the subsonic flow part was especially affected, since as may be seen 

from figure 13 the a line -2.3 still cuts through low values 

of a 0 t especially in the low-subsonic region. Inaccuracies can also 

be noted in the reaching of such higher steady-flow supersonic Mach 
numbers as correspond to values of y = 2 and y = 4. This larger 
error is due to a growth of the error introduced by the rapid 
changes in cross section near the minimum section. Since the error 
increased for increasing values of y, the y = Constant lines for 
higher values of y (up to y = 16.5887) were not plotted in figure 14. 
The line Mq = 5 is, however, included in figure 14 (dashed line) in 

order to indicate the range of the final values for the flow through 
the nozzle designed for a steady-flow Mach number of 5. 

For lower values of y (0.5 and l) the error introduced by rapid 
variation of cross section was still small and the y = Constant lines 
came very close to reaching the required steady-flow Mach numbers. 

(See fig. 14 and table 7.) The gradual approach to the steady-flow 
values (indicated in table V) is the result not only of the approximate 
nature of the calculations but also of the fact that steady flow 
actually presents an asymptotic condition which it would take an infinite 
time to reach. An important conclusion that can be drawn from the 
calculations is that the main changes in flow variables occur very 
quickly and that the remaining changes in the infinite time interval 
are very small. The small number of decimal places used in table V 
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compared to the high number of decimal places used in table TV is 
adequate for the accuracy of the calculations based on the use of 
relatively large steps. Further useful information for the transition 
from unsteady to steady flow may also be obtained by observing the path 
of the lines y = 0.5 and y = 1 in the diagram of a/a 0 against v/a Q . 

The extension of these lines into the area of the final steady-flow ellipse 
signifies that the final steady flow at first tries to establish itself 
at a similar lower speed range, that is, along a similar smaller steady- 
flow ellipse, before it is boosted up to the final steady flow which has 
the critical pressure corresponding to the tank pressure at the mi nimum 
cross section of the nozzle. The boosting process is the reason for 
the consequent zigzagging of the y = Constant lines. Another inter- 
esting effect that occurs in the near— to—maximum speed region concerns 


the speed of the 


v — a 
a 0 


waveB in this region; namely, that they have 


to reduce their speed in the transition from unsteady to steady flow, 
since the unsteady-flow maximum speed is higher (l >a 0 ) than the steady- 

flow maximum speed 


In order to facilitate the interpretation of figure 14, the steady- 
flow Mach numbers corresponding to y values of — O.2325, 0, 0.5, 1 
and 2 were indicated by the M = Constant lines. This auxiliary 
construction is especially useful for y = 0 (minimum cross section) 
since for that case both the initial unsteady flow condition and the 
final steady flow lie bn the M = 1 line. 


Finally, the time required for the transition from unsteady to 
steady flow is discussed. Table Y shows that the values of a<yt 
closest to the steady— flow Mach numbers to be reached lie roughly 
around a 0 t =4.5. From the previous discussion it can be concluded 

that this value will be only very slightly too low in spite of the 
fact that for an inviscid fluid steady flow is reached asymtotically 
after an infinite time. (For the calculations the rounded-off 
value aQt =5 is UBed.) This behavior is in agreement with the fact 

that unsteady disturbances once created in an inviscid fluid do not 
disappear. From a practical viewpoint, though, after a very short 
time the unsteady disturbances become sufficiently small so that they 
are completely dissipated by viscous effects. 


Since the dimension of a Q t is the same as y, that is, a 

length, it can be conveniently expressed in terms of the diameter of 
the minimum cross section of the nozzle. For a nozzle of circ ular 

cross section, both a Q t and y have to be multiplied by \jjr/k 

since d = \jk/ n for F m ^ n = 1. An aQt of 5 then corresponds 

to 5 d\Jit/h . The time it takes to reach steady flow is given 
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Since &q is a constant, the time it takes to reach 


steady flow increases with th 9 nozzle size 
a d m ^ n of 1 foot and room temperature in 

time it takes to reach steady flow will "be 


for a given shape. For 
the tank, for example, the 


~js_ 

1116 ~ 1000 


second. 


The general expression for the transition time of the flow in a 
nozzle with given shape hut varying size may be extended to include 
nozzles of affine shapes (stretched or shortened shapes). This extension 
seems easily understandable inasmuch as the length of a nozzle rather 
than the diameter is the determining factor in th 9 calculation of the 
transition tim 9 , since the length directly affectB the running time of the 
disturbances. The similarity of unsteady— flow phenomena for nozzles with 
the same or affine shape but varying size can also be said to signify 
that a ’’similarity rule” exists for unsteady flows through tubes with 
variable cross section. This fact is of importance to the test engineer. 

Through a process of affine shortening of a given nozzle shape, 
a shape with sudden change in cross section is finally obtained, as 
indicated in figure 15 . For thiB sudden change in cross section, 
steady flow should establish itself after zero transition time. (These 
considerations are based on one— dimensional theory.) On th 9 basis of 
such considerations, the motion of a disturbance of large amplitude 
through variable cross section may be substituted by that through a 
series of cylindric tubes with interspaced abrupt changes in cross 
section. Such a scheme was used in reference 15 without, however, 
giving thorough explanation as to why it should be permissible to use 
the scheme. The nature of the difference between the effects of 
tubes with gradual changes in cross section and sudden changes is 
brought out to some degree in figure 15 . The modifications of the 
original disturbance in this case consist of a weakened disturbance 
continuing to move in the same direction and of a reflected disturbance 
in the opposite direction. Since all the changes of the original 
disturbance take place at the same location y, the continuing weakened 
disturbance as well as the reflected disturbance may be substituted by 
large sudden disturbances (group of dashed centered lines) originating 
at this location. Each centered disturbance may, furthermore, be 
substituted by a single line. From the discussion of the sma l 1 values 
of transition time it can also be surmised that the scheme here discussed 
converges to the characteristic method. 

Finally, a few practical aspects of the short-duration tests are 
discussed. It is important to emphasize that, although practically, 
steady flow conditions will be reached in these very— short-duration 
tests, for most materials commonly used for the construction of the 
test model and the tunnel walls, conditions of temperature equilibrium 
between models, walls, and air flow will not be reached. Because the 
model and the walls will remain close to their initial temperature, the 
amounts of heat transfer through the boundary layer will be different 
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from those encountered in the customary steady-flow tests. The herewith- 
connected changes in boundary— layer behavior (transition and others) have 
to be studied before drawing quantitative conclusions from such tests. It 
should be realized that the constancy of initial temperature of the model 
and the walls considered from a different viewpoint also makes it possible 
to conduct tests with arbitrarily chosen model and wall temperatures at a 
given steady— flow Mach number. Conclusively, it should be added that the 
main difficulty of such very-short-duration tests will be the problem of 
instrumentation. 


PART III.- STUDY OF INVARIANT INTEGRALS OF UNSTEADY FLOW 
DISTURBANCES OF LARGE AND SMALL AMPLITUDES 

GENERAL CONSIDERATIONS 


In part III the behavior of unsteady flow disturbances is treated 
from a viewpoint different from that in parts I and II where a detailed 
study of the complicated pattern of local growths, reflections, and 
re— reflections within the disturbance was made. In part III integrals 
over the whole disturbance invariant with respect to time are developed. 

The value of such invariant integrals is doubtful as long as the 
disturbance, its reflections, and re-reflections belong to the integrand. 
Thus, the main purpose of the present paper is to find integrands which 
permit a separation of the integral effects of the growth of a disturbance 
and its reflections. The logical choice of an integral for the separation 
is one which uses the characteristic parameters of the disturbance as the 
integrand. These characteristic parameters consist of two summands, the 
velocity increment due to the disturbance, and another quantity of the 
dimension of a velocity which represents the pressure increment in the 
disturbance. Depending on whether the sum or the difference of these 
two velocity increments is used, these parameters are associated with 
disturbances traveling with the speed of sound to the right anri to 
the left relative to the fluid. (For more detailed statement see 
section entitled "General Considerations" of part I.) 

The integrals with the characteristic parameters as integrands 
were introduced in reference 5 under the name of pulse areas. In the 
language of the engineer and the physicist, the invariance of this 
integral is also called its conservation. This word causes the question 
to arise immediately as to what relation the conservation of the pulse areas 
would have to one of the three conservation laws, the mass, energy, 
and momentum laws. The relation can be quickly found out from the 
previous definition of the characteristic parameters as a linear 
combination of two velocity increments, one of them representative of 
the pressure increment in the disturbance. It can be easily seen that 
the sum of these two parameters P and Q (or \ and - 4 , respectively) 
will be twice the increment of the flow velocity in the disturbance. 
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Since the integral contains now a velocity increment as integrand, it 
is related to the conservation law for the momentum. It has to he 
emphasized, though, that only a relation exists, hut no identity, 
since the integral of the velocity increment is actually the velocity 
potential. The conservation of the potential here replaces the 
conservation of the momentum which, in contrast to the tube with 
constant cross section, is not conserved because of unknown pressures 
at the inclined walls. Since the invariant integral for the pulse area 
has played such a useful part in expressing an invariance of an unsteady 
disturbance, it is useful to express the invariant integrals for mass 
and energy in terms of the pulse area. These expressions have been 
given in reference 6 without introducing, however, the pulse area and 
applying the restriction of the short disturbance which is necessary 
for the complete separation of the growths (positive and negative) 
and reflections of a disturbance. The reasons for the introduction 
of the short disturbance are given subsequently. 

So far the first step in the process of separating the growth 
in the disturbance and the reflected disturbances lias been shown which 
consists in the introduction of the pulse area. This introduction, 
however, still does not eliminate the re-reflections from the * 
integrand it only labels them and discriminates odd and even numbers 
of reflections. The final step in the separating process consists in 
uncoupling the growth and the reflections of the disturbance and its 
reflections by neglecting the re-reflections during the time a single 
reflection is being produced. This separation (or uncoupling) of the 
original disturbance and its reflection is achieved by restricting 
the length and the amplitude of the original disturbance and the 
inclination of the tube walls to small values (short disturbance). 


ANALYSIS 

Conservation of Potential 


It is shown in references 10 and 14 that for a fluid assumed to 
be inviscid a velocity potential exists for the motion of an unsteady 
flow disturbance of large or small amplitude in a steady flow or in a 
fluid at rest through variable cross section. It is also shown in 
references 10 and 14 that the potential consists of two parts 1 

d0 = v dy - {£- +• iVt (45) 

The next step is to show the invariance of the potential with 
respect to time. This problem is conveniently investigated in the 
coordinate system of y against aot (see fig. 1 6). For the sake 
of discussion a disturbance of large amplitude is assumed to be 
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produced at the time t = 0 in such a manner that behind the disturbance 
its amplitude is zero. The potential 0 is indicated by the area of 
the cross-hatched rectangle. The height of the rectangle presents the 
amplitude of the disturbance, which is, as is shown subsequently, the 
velocity increment through the disturbance. The height, of course, 
does not lie in the y,a Q t plane, but it is conveniently indicated 

in the figure in this plane. The unsteady flow disturbance under 
consideration is assumed to be an isentropic expansion of an arbitrary 
gas produced at t = 0 and to move through an increasing cross section. 

In the course of its motion the amplitude and the length of the 
disturbance change; at the same time reflected disturbances with given 
amplitude and length are given off. In figure 1 6 the reflected potentials 
are indicated by the long, thin cross-hatched areas; whereas the potentials 
of the growing disturbance are shown by short, thick cross-hatched areas. 

In order to show the invariance of the potential 0 with respect to 
time, equation (45) has to be integrated along a path which does not 
enclose the region in which the disturbance is produced, the production 
of the disturbance being an extraneous process (for example, unsteady 
motion of piston, bursting of diaphragm, unsteady motion of side walls, 
combustion or condensation front). The integration along such a path 
is of the following form: 



0 (46) 


Since the disturbance is assumed to move in a steady stream, the 

difference of the integrals with the Bernoulli constant ^io = ^ + i 
as the integrand^ is zero. The remaining integrals are equal to each 
other. The significance of this result is as follows: 


When an unsteady isentropic disturbance of large amplitude travels 
in a steady flow of an arbitrary gas through variable cross section or 
through a gas at rest, the potential of the growing (positively or 
negatively) disturbance itself with the addition of the potential due 
to all its reflections remains constant with respect to time; that 


is, I v dy = Constant. 


The application of this general law to the motion of strong shocks 
through variable cross section seems doubtful, since when entropy 
variations appear in the flow the potential ceases to exist. It would 
seem possible, perhaps, to enforce the existence of a potential even 
for such a case by choosing a path of integration which encloses 
positive and negative gradients in entropy. The situation is analogous 
to that for vortices in steady flow where, at least for skillfully 
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chosen cases (Prandtl’s starting vortex), the proper path of integration 
exists for the maintenance of the concept of potential. 


Although the potential 



has yielded useful conservation laws 


for the "behavior of a large disturbance moving through variable cross 
section, the information obtained from this law is incomplete since 
it only concerns one flow variable and, consequently, does not permit 
a separation of the integral effects of the. disturbances. As indicated 
in the section entitled "General Considerations," the introduction of 

the pulse areas f P dy and f Q dy is a first step toward the 


separation. The parameters P and Q are given for unsteady flow 
of an arbitrary gas through a variable cross section by 


P = v + f(a) 
Q = v - f(a) 


(V7a) 


or for an ideal gas 



(47b) 


Conservation of Pulse Areas 

A conservation law for the pulse areas may be obtained from the 

conservation law 0 = v dy = Constant for the potential 0 by 

expressing the velocity v in terms of the amplitudes P and Q 
given by equations (47). The following conservation law for the pulse 
areas results: - • ■■■ 


(P + Q)dy = Constant 


(48a) 
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or 


P dy + 


Q dy = Constant 


(48b) 


This conservation law for the combination of the pulse areas for the 
disturbance of large amplitude, however, has no direct value for the calcu- 
lation of the behavior of the large disturbance since P and Q (or X 
and p) are interrelated by the system of two nonlinear differential equa- 
tions shown, for example, by equations (6) and ( 7 ) of part I. As stated 
previously, two cases exist for which the integral effects of the 
disturbances can be fully separated. One case is that of short 
disturbances which is discussed in other sections of part III; the other 
is the case of motion of a large-amplitude disturbance in a steady flow 
without Mach number gradient through constant cross section. The 
specification "without Mach number gradient" is made to exclude the case 
of constant cross section at a steady— flow Mach number around one. 

As previously stated, the motion of a large disturbance through 
constant cross section is characterized by the fact that P is 
constant for a disturbance with the amplitude Q, and Q is constant 
for a disturbance with the amplitude P. Thus, the fact is at once 
established that for a large disturbance the individual pulse areas 
are conserved; that is. 


J P d y = Constant (49a) 

Q dy = Constant (49b) 


In order to avoid confusion it is well to point out the difference 
between the symbols P and Q used in reference 5 and the symbols X 
and p used in references 4 and 10. The parameter X and p are 
given for arbitrary gases by 
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or for an ideal gas "by 



The equations are the same as equations ( 4 ) and (5) tut are repeated 
here for convenience. The difference in presentation between P nnH Q 
(equations ( 47 )) and X and n (equations ( 50 ) ) , respectively, maybe 
stated directly in terms of the flow variables a and v for ideal 
gases. Expressed in terms of P and Q (reference 5) 


and 


v = P.. + Q 

2 


(51a) 


. 7 - 1 P-Q 
2 


( 51 b) 


whereas in terms of X and n 


v _ X — u 

^ " 2 • ' (52a) 


and 


_a , 7 — 1 X + ji 

a 0 2 , ~~ 2 ~ 


(52b) 


The convention of- writing the velocity as the difference of tv.', 
flow, parameters X and p.' is taken from steady Uro-d imens i onal 
supersonic flow (see reference 16 ) . The convention of presenting the, 
velocity in terms of the sum of the two flow parameters P and Q is 
taken from reference 17. 
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Finally, a few remarks seem in place concerning the statement made in 
this section that only for the cases of a short disturbance and of constant 
cross section can the integral effects of the disturbances be fully 
separated. It is shown in reference 15 (see, also, discussion in part II) 
that an approximate method of separation exists for a large-amplitude 
disturbance based on substituting the continuous change in cross section 
by a discontinuous change. This approximate method makes skillful use 
of the behavior of a large-amplitude disturbance in constant cross 
section. 


Eelative Growth (Positive and Negative) and Eeflection of the 
Amplitude of a Short Disturbance at a Given Steady— Flow Mach 
Number, Due to a Given Cross-Sectional Gradient or 
Mach Number Gradient 

In the process of determining the growth and the reflection of the 
pulse area of a disturbance the behavior of their amplitudes has to be 
investigated. The equations for the behavior of the amplitudes of a 
small disturbance are given in references 5 and 6. As mentioned in the 
introduction, however, neither of the reference papers f ully deal with 
the growth and the reflection of the pulse area or the amplitude of a 
short pulse. For this reason the derivation in the present paper is 
started with the basic equations for the amplitudes of small disturbances 
moving through a steady flow in variable cross section. The equations 
for air (7 = 1.4) (see reference 5, equations (l4) and ( 15 ) ) are as 
follows : 


§H 


3P + 2Q 


v o + a o + 5 


♦ (i - %)(^} 


lie 

<iy 


(m 2 _ ± )f£ - £ + p ± q + p - o \ 

vb /V 20v o 2M o 


dv c 

4y 


= 0 


( 53 a) 


|s + ( To . a ° . ♦ (1 ♦ ^(^g* 

- (** - - 0 (53(b) 
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For the case of an arbitrary gas the multiples 2 and 3 of P and Q, 

respectively, are given by the expressions - 2-~ . and 

2(7-1) 2(7 -1)' 

respectively, derived from equations (51a) and (51b). The denomi— 

2 _4 


nators 5, 10, and 20 are given by the expressions 
8 


7 - V 7 - V 

and — j For the sake of brevity, equations (53a) and (53b) for air 


are used in the derivation; the quantity 7 is only introduced into the 

final equations. The sums - P - + and — * ^ are equal to the svuns 

Vo + a 0 and v Q — a 0 , respectively. In the case of small disturbances the 
amplitudes P and Q are given by a linear combination of a velocity 
increment and a pressure increment (expressed in the form of an increment 
of the velocity of sound) with respect to steady-flow conditions. 

Since in the expressions for the P and Q pulse (equations (47a) 
and (47b)) the pressure function f(a) has opposite signs, the alternative 
of compression or expansion enters into P and Q with opposite signs. 

The flow velocity (the significant parameter for potential, momentum, 
and pulse areas) enters into P and Q with the same sign. For the 
parameters X and p the situation is reversed. The motion of the 
P pulse is designated in the downstream direction by making the arbitrary 
down-tube direction coincide with the downstream direction. 


Equations (53a) and (53b) are identical with equations (l4) and (lM in 

tfc. equation giveiT for the^ariation 
of the flow variables along characteristic lines for isentropic unsteady 
disturbances of large amplitude given in part I and in references 4 5 

and 10. For the case of small-amplitude disturbances, P and 0 are * 
small and are of the first order. Furthermore, as is indicated by 
equations (53a) and (53b), the first derivatives of P and 0 are 
also small and of the first order. 


For these conditions the following approximations can be made: 

(1) The term may be considered negligible when compared 

with the expression — - + S + P ~ Q 

2Mq 10 ' 

(2) The term 3P j — may be considered negligible when compared 
with v Q + a Q . 

(3) The term — * 3Q may be considered negligible when compared 
with v Q — a Q , with certain qualifications. 


The qualifications mentioned in approximation (3) 
fact that as approaches one, v Q - a Q approaches 


are due to the 
zero. If v Q - 


o-» 
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however, approaches zero, the term - jr 


may not he neglected. In 


view of these conditions, the behavior of the Q pulse and its reflections 
has to he investigated separately for the case of 1$, = 1. The assurance 

whether or not the term —22 be actually neglected near Mq = 1 

may he obtained by comparing the results of the separate investigation 
at = 1 with the results obtained by neglecting — + _ 2 ., 


The approximations made so far in equations ( 53 a) and (53b) are all 
a direct consequence of the use of disturbances of small amplitude. The 
approximations, however, still do not simplify the simultaneous differ- 
ential equations (53a) and (53b) sufficiently to make the solution easily 
interpretable. For a simple solution of the simultaneous equations, a 
type of uncoupling of P and Q would be desirable. More specifically, 
the uncoupling would signify that in the equations ( 53 a) and ( 53 b) P and Q 
could be alternately neglected depending on whether one is attempting to 
solve for the growths of the pulse or for the reflected pulses. The 
"uncoupling" process was substantially used in reference 5 without explaining 
its full meaning. For a clear understanding of all approximations made in 
the calculations, the concept of uncoupling is explained briefly: For 

example, as a Q pulse moves through a tube with variable cross section, 
its amplitude Q grows to Q ± dQ and gives off a reflected pulse ±dP. 

The building-up process of the reflected pulse P (traveling with the 
speed v 0 + a Q ), during its motion through the growing Q pulse to its 
full strength dP, starts from zero immediately behind the head of the 
Q pulse and ends when the reflected pulse at its full strength dP 
leaves the pulse Q. Because of the building-up process, P is at first 
very small compared with the average Q. The smallness of P compared 
with Q can be kept forever when the Q pulse has a length dy 

and when the gradient of cross section dF/dy or, consequently, 

dJF 

^ = dy dyj the change in cross section of the tube occupied by the 

disturbance Q, is small. A disturbance for which amplitude, length, and 
cross-sectional gradient are subject to these restrictions is called a 
, short disturbance (or pulse) in the present paper. For the same reason, 
for the growing of the pulses to . Q ± dQ, P may be neglected compared 
with Q. The analogous considerations are true for P ± dP and the 
reflections dQ, where Q may be neglected compared with P. 

Before attempting to solve equations. ( 53 a) and (53b) for the 
desired quantities, the equations are written in different form by 
performing the following transformations: 



+ a. 


3 P + 2 Q 


!>*■(• 


+ a. 


3P + 


2 Q|dP 

<iy 


( 5 ^a) 
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and 


I * (v„ - - (vo - *0 ♦ (5M0 


At first equations (54a) and (541)) are substituted in equations (53a) 
and (53^); use is made of approximation (3) without qualifying 

assumptions for Mq = 1, that is, with neglect of 3Q compared 

with v Q - a Q . 


Since presenting the deformation and the reflection of pulses in 
terms of the variation in cross section is desirable, the following 
substitution obtained from the continuity equation for steady flow is 
made additionally: 


dv o v o 1 dJT 

dy % 2 - 1 F dy 


(55) 


The substitutions transform equations (53a) and (53b) into 


( v ® + a o)^ + 


v o 1 dF /, 
^2 _ i F dy [V- 1 - 



K 2 - 1) 


(l ± Q 

V 2% 


+ 



= 0 


(56a) 


and 


( v o 



v o 1 d LF 
- 1 F dy 





P + 0, 

2Mo 



(56b) 


The fact that dy may be canceled in equations (56a) and (56b) has 
the significance that the equations are independent of the scale of y. 

In other words, the equations for short disturbances may also be 
interpreted to apply to small— amplitude disturbances of arbitrary length 
traveling through a small discontinuous change in cross section. As for 
the case of short disturbances ho interference exists while the reflection 
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is "being produced; the lack of interference is due to the fact that the 
reflection is produced in the zero length of the discontinuity and that 
before and behind the small discontinuity the cross section is constant. 
This last interpretation of the approximation of short disturbances 
establishes a connection with the approach given in reference 15. The 
behavior of a small-amplitude disturbance of arbitrary length traveling 
through a small continuous variation in cross section may be obtained by 
pairing small discontinuous changes together in such a mann er that the 
distances between them equal the lengths of the growth or the reflection 
of the small disturbances. The variation of the lengths of disturbances 
with Mach number is discussed in a subsequent section. 


The following relations can now be deter mine d - 

dP/P represents the relative growth dP of the pulse P. 

dQ/P represents the relative reflection dQ produced by the 

pulse P. 

dQ/Q represents the relative growth dQ of the pulse Q. 

dP/Q represents the relative reflection dP produced by the 

pulse Q. 


The following expressions for an -arbitrary gas are obtained from 
equations (56a) and (56b) by making use of equations (4o) to (42.) of 
part II and the previously developed "uncoupling" effects, that is, 
alternately neglecting P and Q: 


dP _ _ 1 1 + 7 g Mo 2 dP 

P ~ 2 (Ho + l) 2 P 

1 - 1 ~ 1 M 2 
dQ _ _1 _ 2 ^ dP 

P " 2 (Mo -l ) 2 » 


dQ _ _ 1 1 + 7 2 ^ dP 

Q ' 2 (Mo - l ) 2 P 


( 57 a) 


(57b) 


( 57 c) 


dP _ _ 1 1 ~ 1 Ha 2 dp 

Q " 2 (Mq + l) 2 P 


( 574 ) 
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Equations (57) reduce to the acoustic case for Mg = 0. Expressing 
equations (57) also in terms of the steady-flow Mach number gradient is 
useful for the remaining investigations in this paper. This is done by 
the substitution of 


dE _ 1 - Mp 2 dfr^ 

F ' 1 + 2-f-i *> 

obtained from the steady— flow continuity equation into equations ( 57 ). 
The following equations result 


dP 1 - dMp 

P 1 + Mg 2Mg 


(58a) 


y _ 

dQ = Mp + 1 1 7 T- tyj 2 dMg 

P ‘ ^0 ~ 1 1 + 2Lzl 1 ^2 2Mg 


(58b) 


dQ _ % + 1 dMo 

Q Mg - 1 2Mg 

, i y 2 

dP = 1 - Mp Z 2 to dMg 

Q " 1 + ^ 1 + Z_I_i Mp2 2Mo 


(58c) 


(58d) 


A closer study of equations (57) and (58) reveals first of all the 
interesting fact that the amplitudes of the relative reflections 
expressed by the ratios dP/Q and dQ/P change their signs at the 

supersonic steady-flow Mach number Mg = which equals \{J for 

air; that is, at Mo = \jy ^ both the P and the Q pulses give off 

no reflections. A more detailed study of equations (57) and (58) is 
given subsequently. 

Equation ( 57 d) indicates that, for example, for the upstre am travel 
of a Q pulse in the subsonic part of a diffuser (0 < Mg < l) the 
reflected pulse dP is negative since it moves in the downstream 

direction 0p > (A from zero (while it is being produced from zero). 



56 


NACA TN No. 1878 


Since the factor 1 - - - 1 > 0 for 0 < ^ < 1, the sign of dP 

consisting of the product of one negative quantity and two positive 
quantities, that is, the product (-)(+)(+) is negative. In order to 
determine whether the reflection dP is of the expansion or compression 
type if Q is, for example, a compression, equations (Uj) must he 
considered. - They state that P and Q are of the opposite sign. Thus, 
if according to equation ( 57 d) dP has a sign opposite to Q, the 
additional information of equations (^7) indicates that dP and Q 
are of the same type; that is, if Q is a compression, dP is a 
compression also. For Mq = 0, this result is in agreement with the 
well-known acoustic behavior. The result obtained may also be conveniently 
discussed in terms of equation ( 58 d) which uses the Mach number gradient. 
For 0 < Mq < 1, the product of signs (+)(+)(—) is negative (dM 0 < 0 ) 
which naturally agrees with the result obtained from equation (57d). 

Another especially interesting illustration of equations (57) and (58) 
concerns the reflections dQ produced by a P pulse traveling down- 
stream in the supersonic part of a nozzle. In this case, both the 
pulse and its reflection travel down the tube since for supersonic 
flow an upstream motion with the speed v 0 - a c is a down-tube motion. 

The behavior of the reflections of a P pulse is described by 

equations ( 57 b) and (58b). For 1 < Mq < 1 the sign 

product (—)(+)(+) is negative for equation (57b) and (—)(+)(+) is 
negative for equation ( 58b ) . Since P and Q are of the opposite 
sign, both t he original pulse and the reflection are of the Same type. 

For < Mq < » the sign products, for ( 57 b) and (58b) are 

both (-)(-)(+) which is positive; thus, in this case the amplitude of 
the original pulse and that of the reflection are of opposite types. 

With the aid of the sign products the types of the growths wnfl 
reflections of the downstream-moving P pulses and the upstream- 
moving Q pulses not mentioned so far may be also determined from 
equations ( 57 ) and (58) with two exceptions. The two exceptions concern 
the behavior of the reflection dQ produced by P (equations ( 57 b) 
and (58b)) and the growth dQ produced by Q (equations ( 57 c) and (58c)) 
at a steady-flow Mach number of Mq = 1 . At that Mach number, dQ becomes 
infinite with respect to both the amplitudes of the pulse P and the 
pulse Q. This growth toward infinity is in contradiction to the 
assumption of small disturbances for- which this result is obtained. 
Equations ( 57 ) and (58) only hold for a range in which the order of 
magnitude of "small" is not exceeded for the amplitude dQ. The range 
can be increased such as to include the closest proximity of =1 
by choosing values close to zero for the amplitudes P and Q of the 
pulses. The deterioration of dQ at M = 1 itself, however, cannot 
be eliminated. Similar considerations apply to the gradient bQ/by. 
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Another case concerning the behavior of the amplitude of a small 
disturbance near a steady— flow Mach number of one concerns the steady- 
flow equilibrium conditions involving a shock near = 1 . A 

discussion of the equilibrium of shocks, however, is outside the scope 
of this paper. 

The discussion so far gives a. very unsatisfactory picture about 
the singular behavior of a small disturbance at ^ .= 1. This picture 
does not seem to show anything significant, for it does not seem 
reasonable that an important quantity like the energy or the mass of 
a reflected disturbance dQ produced by a downstream disturbance P 
should become infinite at 1^=1. The reason for the discrepancy 
between the physical expectations and the results obtained so far 
lies in the fact that only amplitudes of the small disturbances have 
been discussed so far and not the integrals over the whole disturbance 
extension as required for the mass, energy, and pulse area. A 
discussion of these integral quantities is given in the following 
section. 

Before going into these problems, the behavior of the amplitude dQ 

near Mq = 1 is explained, without neglect of the term — + ^ 

/ » 5 

compared with v Q — a D in equation ( 53 b). For a short disturbance 
at = 1, equation (53b) yields: 


dQ 

Q 



if P is neglected (uncoupled) with respect to Q. The neglect 
(uncoupling) of Q with respect to P (equation ( 53 b)) yields 


( 59 ) 



(60) 


By comparing equations (59) with equations (57c) and (58c) and 
equation (60) with equations (57b) and (58b), it may be seen that 
similarly to equations (57c), (58c), (57b), and (58b), equations (59) 
and (60) have an infinity. The only difference is that the infinity 
of equations (59) and (60) does not occur at = 1, but rather at a 

Mach number corresponding to a velocity v Q — a G + Similar 

considerations can be made for equation (6l) and equations (57b) 
and (58b). The physical meaning of . these more exact equations is 
that the growths or the reflections, of the Q. pulses do not actually 
accumulate exactly at the minimum cross section of the nozzle. 
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Conservation of Mibb, Energy, and Sum of Pulse Areas 
for Short Disturbances 

The conservation of mass and energy for large and small disturbances 
traveling through tubes with variable cross Bection without the benefit 
of uncoupling between the original disturbance and its reflection is 
discussed in reference 6 . For the case of short (uncoupled) disturbances 
it is clear that the amount of mass, energy, and pulse area of the growth 
of the original disturbance equals the amount of mass, energy, and pulse 
area, respectively, traveling in the reflected disturbance, since no 
mass, energy, or pulBe area can disappear by the interference between 
the original and the reflected disturbance during the building-up process. 
The proof of the conservation for short disturbances is given subsequently 
for the simplest case of the pulse area. For the sake of convenience, 
equation (48b) is rewritten here: 


P dy + 


r 

/ Q dy = Constant 


Under the assumption that the growth and the reflection of the pulse 
area have the shape of a rectangle with the lengths Ayp and A yq 

equation (48b) can be presented in the following form: 


P Ayp + Q A 7 q = Constant 


( 6 la) 


or 



+ 



(6lb) 


The differentiation in equation ( 6 lb) results in the following 
equation: 


dP Ay p + P <3Ay p + dQ Ayq + Q d/y-Q = 0 


( 62 ) 


If one chooses, for example, to express the conservation law of 
a relative growth dP/P and the relative reflection dQ/P, for the 
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case of a short P pulse, Q starts from zero and the last term can 
he neglected. The conservation law thus assumes the following form: 


dP £y p + P dAy p = -AVqdQ (63a) 


or 


d(P 47 P ) = -Ay Q dQ (63b) 

Equation (63b) indicates that the pulse area of the growth d(P Ay p ) 
equals the reflected pulse area, since P and Q are of opposite Bign. 
Equation (63a) may also be expressed in terms of th9 relative amplitudes 
given in equations (57) and (58): 


dP dAyp _ dQ 

P Ay p Ay p P 


( 64 ) 


Pulse Length and Pulse Time for Short Disturbances 

It was stated previously that the travel of a short disturbance 
through a small continuous change in cross section is identical with 
the travel of a small disturbance through a small discontinuous c hang e 
in cross section. For the case of a small discontinuous change in 
cross section (see fig. 17), the pulse times Atp and AtQ are equal 


and the Mach number gradients are equal and opposite, 
are expressed by the following equations: 

These conditions 

Ay p = At p (v 0 + a Q ) 

( 65 ) 

Ay Q = At Q (v G - a Q ) 

(66) 


It should be recalled that without the restriction to a short 
disturbance, reflections and re-reflections occur which will necessitate 
a step-by-step integration for the determination of the pulse lengths 
and pulse times. 
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Relative Growth and Reflection of the Pulse Area of a Short 

Disturbance Due to a Given Mach Number Gradient, as 

a Function of the Steady— Flow Mach Number 

The behavior of the relative reflection d(Q Ay Q )/P Ay p of the 
pulse area P £y p and of the relative growth d(Q A/q)/Q of the 

pulse area Q Ayq is investigated first since their amplitudes were 

shown to degenerate for Mq = 1. The relative reflection may be written 
as follows* 


d ( Q Ago) _ dQ Ay Q Q dAjrq 

P p Wp + T 5 Ay p 


(67) 


The second term is small as Q starts at zero; therefore, 

d (Q 4 r Q ) dQ A[Q 

P Ay p P Ay p 


( 68 ) 


The terms dQ/P and Ayq //\y P are given in terms of the Mach number. 
The relative reflection dQ/P is expressed already in equation (58b). 
The ratio Ayq/Ayp is formed by substituting equations ( 42 ) and ( 4 l) 
of part II into equations (65) and (66), respectively, which gives 


£y p = aoAt p 


Mq + 1 


and 


(69a) 


= a oA t Q 


- 1 


A + ^ ^ 


(69b) 


The required ratio of Ay 0 and Ay p is then given by 

AyQ = v o ~ a o Mq - 1 
Ayp v o + a o m 0 + 1 


(70) 
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0 

The value of the relative reflection is then 


d (Q Arp) s H, . 1 1 ~ ^ ^-1 dH, 

p Arp Mo - 1 1 + 2 L =-1 n,2 H, + 1 2^, 

1 ~ Ar 2 % a dH, 

1 + ^ vtf 2M ° 


Thus, the magnitude (not the shape) of the relative reflection of the 
pulse has no singularity at Mq = 1. The relative growth d(Q 4 v q )/Q Ay 0 
may he written as follows: ' ^ 


d(Q Ayg) _ d Q dAyg 

Q 4 yQ Q + AyQ (72) 


The quantity dQ/Q in terms of Mach number is already expressed in 
equation (58c). The ratio dAyq/Ayq is obtained by logarithmic 
differentiation of equation (69b) 


dAyQ 

47 q 


log Ay q = 


dM o 

Mo - 1 


7-1 


f^d^ 


1 + V Mo 2 


The value of the relative growth is then 


d ( Q ^q) _ Mq + 1 dMo df^, ^ dM o 

Q Ay Q Mq - 1 2Mo + Mo - 1 ~ ~ y — 1 p 

1 + — ^ — Mo" 


(73) 


L 2 Jl 


( 7 M 


Thus the relative growth of the pulse area has no singularity at Mq = 1 . 
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Next , the "behavior of the relative reflection d(P Ayp ) fq AyQ and 
of the relative growth d(P Ayp)^P Ayp of the pulse area is investigated. 
The relative reflection may "be written as follows: 


d(p Ayp) ^ ^p Ayp p dAyp 
Q A7"q ~ Q Ay-Q + Q A7 q 


( 75 ) 


For short disturbances (P neglected) the relative reflection is 

d(P Ay p ) dp Ay P 
Q Ay"Q “ Q Ay Q 

Equation (76) is written in terms of Mach number by substituting 
equations (58d) and (70) 


d(P Ay P ) _ 1 - ^ 1 ~ ^~2~ Mo 2 % + 1 d% 
Q A7Q 1 + p + 7-1 ^2 Mo - 1 2 Mq 


^ 2 ^° 2 dMj 

! + 2^1 *2 2* 


( 77 ) 


The interesting fact is indicated by equation (77) that the relative 


reflection of the pulse area 


Q is zero at 



which 


equals \[J for air but will not be zero at Mo = 1, in spite of the 

fact that the amplitude is zero there. In reference 5 it is concluded 
from the zero amplitude that the relative reflection of the pulse area 
would have to be zero at Mo = 1 . In the next section of the present 
paper a numerical evaluation of this inaccuracy is given. The relative 
growth d(P ATp)^P Ayp be written as follows: 


d(P Ay p ) _ dP dAy p 

P Ayp P Ayp 


( 78 ) 


Equation (78) is expressed in terms of Mach number by substituting 
equation (58a) and the logarithmically differentiated form of 
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equation (69a). The logarithmic differentiation of equation (69a) 
results in 


Afp 


d 


log Arp = 


1 + Mo 


7-1 




1 + 2L vr J: Mo 2 


The relative growth then is 

d(P Af P ) 1 - d^ 

PAF P = 1 + I^2^, + 1 + ^ ~ ~ 7 - 1 


d Mo 1U t k 

2 Mo 1 + ^-1 ^2 


(79) 


( 80 ) 


Thus , the relative growths of the pulse areas P AFp (equation (80)) 
and Q AFQ (equation (7b)) as well as the relative reflections 
(equations (77) and (71)) are th9 same. 


Accumulated Growths and Reflections of the Pulse Area of a Short 
Disturbance as a Function of Steady— Flow Mach Number 

The accumulated growths and reflections are obtained by integrating 
the relative growths and reflections, respectively. The integration of 
the relative growths d(P AFp)/P AFp and d(Q AFq)/Q AFq can be 

performed without difficulties. The integration (accumulation) of the 
relative reflections, however, cannot be presented in closed form for 
the following reasons: The short disturbance during its motion through 

variable cross section gives off reflections, the accumulated value of 
which will soon have a length larger than that required for a short 
disturbance; furthermore, within this accumulated reflection, re-reflections 
wil*. occur. Consequently, the integration of the reflections would 
require the use of a point-by-point method like that given in part I. The 
fact that the values of the accumulated reflections cannot be represented 
by an integration in closed form 1 b not as unfortunate in the present 
case as might be expected since a conservation law exists for the pulse 
area which permits determination of the acc umula ted relative reflection 
readily from the easily integrable relative growth. In the following 
discussion the variation with Mach number of the accumulated relative 
growths is given. The integration of the relative growths d(Q AFq)/Q AFq 

and d(P Ayp)/P AFp, as given in equations (74) and (80), results in 
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or 


i (Q too) 

QA5^q 



d(p Ayp) 

P Ar P 

► 


log^Q Afq) = log(P A7p) 



(81) 


The definite integration from a reference pulse area 

to Q ATq and the substitution of the definite integral J Q OLy over the 
pulse length ATq at a given time t for Q ATq result in 



Q dy = Cq 



(82) 


The definite integration from a reference pulse area Pq^p^ 

to P Z^Tp and the substitution of the definite integral J* T dy over the 
pulse length /ty-p at a given time t for P Z$rp result in 



Equations (82) and (83) are formally equivalent, the difference in 
constants of the two expressions is due to the fact that the P pulse 
moves with the velocity v Q + a 0 , whereas the Q pulse moves with the 

velocity v 0 — Uq. In terms of the Mach numbers appearing in 
equations (82) and (83) the meaning of this difference is that the 
Mach number in one equation will have the negative value of that in the 
other. The constants in equations (82) and (83) are thus related by 
the equation 

C 2 = «J C 1 


The significance of the imaginary factor J in the equation is that 
a real P pulse can never become a real Q pulse. 

It. may be seen from equations (82) mid (83) that the accumulated 
relative growths of the pulse areas Q A^q ' and P Z\y p have a maximum 

for Wo = ]fj~i which equals {5 for air. The accumulated reflections 

thus remain constant only for Mq = which equals \f$ for air. 
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(The relative reflection is zero.) Equations (82) and (83) permit a 
numerical evaluation of the variation of the accumulated growth with 
steady— flow Mach number (see table Vi) . Table VI indicates that the 
assumption made in reference 5 concerning the conservation of pulse 
area at = 1 is still fairly accurate ^ .since f(^>) representing 

the variation with Mg of the accumulated growth remains fairly constant 
in the range of M Q = 1 to = 5 . The main changes in f(Mo) occur in 
the regions of low subsonic and high supersonic Mach numbers. 


Accumulated Growths and Reflections of the Mass of a Short 
Disturbance as a Function of Steady— Flow Mach Number 

For reasons similar to those in the case of the pulse area only 
the accumulated growths are calculated. The basis for the calculation 
is the following definite integral over the pulse length Ajr at a given 
time t (see also reference 6): 


I 


pF dy 


m 


The quantity p represents the excess in density compared with steady 
flow at a given cross section F. The first step in the separation 
of the growth and the reflection is obtained by expressing p in terms 
of P and Q (see equations (bj)) 


P = v + — - — a 
7-1 


Q = 



or 


P + Q = 2v 


P - Q = — a 
7-1 

The next step in the separation is obtained by the introduction of 
the assumption of short disturbances. Since for the pulse area the 
variations of accumulated growths of pulses P and Q gave equations 
of the same form (equations (82) and (83)), the behavior of one of the 
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pulses is only treated here. For the growth of, for example, the 
pulse P, Q is neglected; that is. 


= 2v = 


7-1 


( 85 ) 


The quantity a is expressed in terms of p for a htdaI 1 disturbance 
by the following isentropic relation: 


_P_ = 2. _a 

P 0 7 - 1 a 0 

The substitution of equation (86) into equation (85) results in 


( 86 ) 



Substituting equation (87) into equation ( 84 ) gives 


( 87 ) 


P pF dy 


u • 



PQ^O 

2a 0 v 0 



( 88 ) 


where PqFv 0 /2 has a constant value, let us say, C 3 according to the 
continuity equation for steady flow. The equation (88) is written in 
terms of Mach number by substituting for v 0 and a Q the values from 


equations ( 40 ) and (39) in part II and for 



the expression from 


equation (83). The following equation results: 



C3C2 




( 89 ) 


The function 



is tabulated in table VII. 


At M 0 = \JyzTi which equals \p5 for air, f (M q) has a minimum. The 

physical significance of this behavior is that the accumulated growth 
of the mass of a short disturbance and, consequently, the reflection 
will be zero at that Mach number. No special behavior occurs for Mq = 1. 

The infinity indicated in table VII for 1 ^ = 0 (acoustic case) is due 
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to the fact that from the point of view of the steady— flow problem, 
the acouBtic case has infinite cross section. The infinity at ^ = » 
is due to the same reason. 


Accumulated Growth and Deflections of the Energy of a Short 

Disturbance as a Function of Steady— Flow Mach Number 

The basis for the calculation of the accumulated growths is the 
following definite integral over the pulse length £y at a given 
time t: 


1 

1 


( E o + E )K + e) E 


dy - Je q 


dy 


which when the second— order terms are negligible compared with first- 
order terms is equal to 


/f P o E dy + J: 


FE_p dy 


(90) 


In the above expressions, E + E 0 is the total convective energy 
stored per unit mass of the small disturbance 


E + 


•nQ 

E o = V + c v T - e Q 


= 7^ A 2 

2 7(7 - 1) 


e 0 


(91) 


and E Q is the contribution of the steady flow to the convective energy 
of the sma l l disturbance 


E o = ^§“ + c v T o “ e 0 (92a) 

or, since E Q is a steady— flow contribution, it may also be written in 
the form 



(92b) 


The quantity eo is an arbitrary constant, which enters equation (91) 
because the convective energy is a potential and, therefore, can be 
referred to an arbitrary level. The terms V, T, and A are the total 
values of the velocity, absolute temperature, and velocity of sound. 
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respectively. Neglecting the second-order terms, which are negligible when 
compared with the first— order terms, in equation (91) gives 


2aoa 


E = v 0 V + . 

° 7(7 - 1) 


(93) 


The substitution of equation (85), (87), (92b), (93)* (39)* 
and (40) into equation ( 90) results in 


1 + ^2 


Fp 0 E dy + J FE 0 p dy = C3 / P dy + — '' c 3 I p d y 


Co 1 + 

-e. —2 E / p dy 

0 a 0 2 17 


(9<0 


The substitution of equation (83) into equation (9^) gives 


Fp ° E 47 + J ra ° p dy - C & U , ^ 1 , j + r±4 


!1 + Mp 2 




e 0 \ 
ao 2 ' 


7 - 1 „ 2 
“2— % 
% 


(95) 


For reasonable values of the constant 6 q, the accumulated growth 
of the energy of a short disturbance shows a behavior similar to that of 
the pulse area; namely, it has an extremum at M 0 = ^ 2 ^ which 
equals \[5 for air and shows no special behavior at Mo = 1. 


CONCLUSIONS 


Part I .— The point-by— point method developed for the calculation of 
unsteady flows through tubes with variable cross section permits a simple 
presentation of* the interaction of strong shocks and large temperature 
contact, discontinuities, a detailed treatment of which had not been given 
so far. The point— by— point method permits furthermore a presentation of 
shock calculations and calculations of flows with initial entropy gradient 
in a form convenient for computation by use of computing machines. Under 
certain assumptions the formulas established inay also be used for the 
calculation of flows with continuous heat addition over a large space. 
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Part II .- The calculations, made of the flow pattern created "by 
bursting into a vacuum of a diaphragm at the minimum section of a 
supersonic nozzle without a second' throat, indicate the following: 

The transition time from .the starting of the flow to the attainment 
of approximately steady flow conditions is sufficiently short to 
permit the use of very-short-duration tests. The transition time for 
the specific nozzle is presented in such a form that a "s imi larity 
rule" can he established concerning the transition time for nozzles 
of different size but of the same or affine shape. 

Part III .- The restriction to short disturbances permits a 
simple presentation in terms of steady— flow Mach number of the growths 
and reflection of pulse area, mass, and energy of a disturbance 
traveling through a steady flow in a tube with variable cross section. 
The calculations show the interesting result that the conditions for 
zero reflection of mass, energy, and pulse area exist at a steady-flow 

Mach number Mq = \ j^ j ^ (which equals 0 for air), where 7 is the 

ratio of specific heats, rather than at 1 ^= 1 ; they also show that 
for practical purposes for the range of Mach number from 1 to 5 the 
reflections are small enough so that the mass, energy, and pulse area 
of the original disturbance may be considered constant. However, at 
low subsonic and at high supersonic Mach numbers, the reflections may 
not be neglected. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va. , January 14, 19^9 





The differentiation results in 






The analytic expression for n, obtained by differentiating equation (Bl), is 
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TABLE I 


V AEIATIOW OF Bi^j WITH Av j f An Jor^y 


AND IT FOR A SHOCK IN AIR 


[7 = 1.4] 


Av/aj_ 

ag/a-L 

Au/a^ 

n 

A v / a l 

®2/ a l 

Au ja^ 

IT 

0 

.0330 

.0654 

.0972 

.1284 

•1591 

.1893 

.2190 

.2483 

.2771 

.3056 

• 3336 
.3613 
.3886 
.4156 
.4423 
.4687 
.4948 
.5206 
•5461 

• 5714 
.5965 
.6213 
.6459 
.6703 

1.0000 

1.0066 

1.0131 

1.0195 

1.0257 

1.0320 

1.0381 

1.0442 

1.0502 

1.0561 

1.0621 

1.0680 

1.0738 

1.0797 

1.0855 

1.0913 

1.0971 

1.1028 

1.1086 

1.1144 

1.1201 

1-1259 

1-1317 

1.1374 

1.1432 

1.00 

1.02 

1.04 

1.06 

1.08 

1-10 

1.12 

1.14 

1.16 

1.18 

1.20 

1.22 

1.24 

1.26 

1.28 

1.30 

1.32 

1.34 

1.36 

1.38 

1.4o 

1.42 

1.44 

1.46 

1.48 

1.00000 

.99990 

.99990 

.99970 

•99950 

.99890 

.99830 

.99721 

.99602 

.99453 

.99285 

.99069 

.98834 

.98571 

.98270 

-97934 

.97570 

.97182 

.96759 

.96302 

.95822 

.95302 

.94760 

.94197 

.93598 

0.6944 

.7184 

.7422 

.7658 

.7892 

.8125 

.8356 

.8585 

.8813 

.9040 

.9265 

.9488 

•9711 

•9932 

I.0152 

1.0370 

1.0588 

1.0804 

1.1020 

1.1234 

1.1447 

I.1660 

1.1871 

1.2082 

I.2291 

I.2500 

1.1490 

1.1548 

I.1606 

1.1664 

1.1723 

1.1781 

1.181a) 

1.1899 

1.1958 

1.2017 

1.2076 

1.2136 

1.2195 

1.2255 

1.2315 

1.2376 

! 1.2436 
1.2497 
1.2558 
1.2619 
1.2680 
1.2742 
1.2804 
1.2866 
1.2928 
1.2990 

1.50 

1-52 

1.54 

1.56 

1.58 

1.60 

1.62 

1.64 

1.66 

1.68 

I.70 

1.72 

1.74 

1.76 

1.78 

1.80 

1.82 

1.84 

1.86 

1.88 

I.90 

I.92 

1.94 

1.96 

1.98 

2.00 

» ©<3 ssssi M a Kite Set 




m 


Ay / a-i 

d(a 2 /a 1 ) 

dfAv/ajL) 

0 

0.2000 

.05 

.2003 

.10 

.2011 

.15 

.2024 

.20 

.2041 

.25 

.2062 

.30 

.2089 

.35 

.2120 

.40 

.2152 

.45 

.2189 

.50 

.2229 

.55 

.2270 

.60 

.2316 

.65 

.2363 

.70 

.2412 

.75 

.2462 

.80 

. 2514 : 

.85 

.2567 

.90 

.2620 

.95 

• .2674 

1.00 

.2729 

1.05 

.2839 

1.10 

.2890 

1.15 

.2948 

1.20 

.3008 

1.25 

.3073 

8.00 

.5106 
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TABLE III 


VARIATION OF — — = n WITH &2/& ± 

d^/aj 

FOR A SHOCK IN AIR 


0 

-.0041 

-.0176 

-.0380 

-.0652 

-.0980 

-.1367 

-.1786 

-.2243 

-.2724 

-.3254 

-.3777 

-.4316 

-.4883 

-.5430 

-.5982 

-.6537 

-.7079 

-.76X9 

-.8145 

-.8649 

-.9151 

-.9625 

I.OO81 

1.0522 


1.1490 

1.1548 

1.1606 

1.1664 

1.1723 

1.1781 

1.1840 

1.1899 

1.1958 

1.2017 

1.2076 

1.2136 

1.2195 

1.2255 

1.2315 

1.2376 

1.2436 

1.2497 

1.2558 

1.2619 

1.2680 

1.2742 

1.2804 

1.2866 

1.2928 

1.2990 


d(a 2 / a i) 

- 1.0944 

- 1.1347 

- 1.1720 

-1.2076 

- 1.2413 

-1.2720 

- 1.3015 

- 1.3286 

-1.3528 

- 1.3756 

- 1.3959 

- 1.4145 

- 1.4311 

-1.4458 

-1.4583 

- 1.4696 

-1.4785 

- 1.4862 

- 1.4923 

-1.4967 

- 1.4993 

-I.5OO9 

-1.5012 

-1.5000 

- 1.4977 

- 1.4940 
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TABLE 71 


VARIATION WITH MACH NUMBER OF THE FUNCTION 


f(M 0 ) = 




Mo 2 


REPRESENTATIVE OF THE ACCUMULATIVE GROWTH 
OF FUISE AREAS FOR AIR 


[7 = i.t] 


M 

0 

f(Ho) 

0 

0 

.1 

.316 

.3 

.543 

• 5 

.690 

• 7 

.798 

1.0 

.913 

1.4 

1.002 

1.5 

1.017 

1.6 

1.029 

1.8 

1.045 

2.0 

1.054 

\I5~ 

1.057 

2.4 

1.056 

2.5 

1.054 

3.0 

1.035 

4.0 

.976 

5.0 

.913 

10.0 

.690 

100.0 

.224 

00 

0 
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TABLE HI 


VARIATION WITH MACH NUMBER OF THE FUNCTION 



REPRESENTATIVE OF THE ACCUMULATIVE GROWTH OF MASS FOR AIR 

[ 7 = 1 . 4 ] 


M o 

f(M c ) 

’ 0 

00 

.1 

3.165 

• 3 

1.8142 

•5 

1.449 

•7 

1.252 

1.0 

1.095 

1.4 

•997 

1.5 

.983 

1.6 

•972 

1.8 

•957 

2.0 

.949 

{ 5 

.946 

2.4 

.947 

2.5 

.949 

3.0 

.966 

4.0 

I.025 

5.0 

1.095 

10 .0 

1.449 

100.0 

4.473 

00 

0 ° 

1 
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Figure 6.- 


Variation of 


drr 

d ( q.o /a-i^ 


n with a 2 /a 1 for a shock 
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Figure 7.- Sample results of bursting of a diaphragm presented in a 
diagram of y against a^t. 
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Figure 12 . Presentation of disturbance velocities in the plot of \[ 2 i 

against v. 









ing through sudden change in cross 


e plot of y against a^t. 
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Figure 17.- Motion of a small disturbance through a small discontinuous 

change in cross section. 




